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soft  silicon  rubber  cable. 

The  first  experiment  used  cables  6  feet  long.  Five  dis¬ 
tinct  geometries  were  measured,  including  simulation  of  the 
anchor-^ast  deployment  of  a  single  mooring  leg,  the  relaxation 
of  a  subsurface  mooring  from  a  displaced  condition,  and  the  re¬ 
sponse  of  a  load  suspended  along  a  cable  fixed  at  one  end  and 
moved  around  a  circle  at  the  other.  The  deployment  and  relaxa¬ 
tion  simulations  were  repeated  in  the  second  experimental  series 
using  cables  60  feet  long.  The  tension  was  measured  during  the 
third  series  while  a  weight  suspended  from  the  nylon  or  rubber 
cable  was  paid  out  or  reeled  in  from  a  small  winch.  The  winch 
base  could  be  oscillated  vertically  to  simulate  wave  action. 

'The  fourth  experiment  was  the  deployment  of  a  full-size 
instrumented  subsurface  mooring  in  2,500  feet  of  water  off  Kauai 
Hawaii  as  part  of  the  joint  Mooring  Dynamics  Experiment  (MDE) 
conducted  in  late  1976.  * 

SEADYN  is  a  general  three-dimensional  model  of  the  dynamic 
response  of  cable  networks  to  environmental  changes  using  a  non¬ 
linear  finite  element  techniques  in  the  time  domain.  SNAPLOAD 
uses  the  lumped  parameter  method  to  model  the  dynamics  of  seri¬ 
ally  connected  cables  suspended  in  a  vertical  plane.  Both 
models  are  shown  to  reproduce  all  the  significant  motions  and 
^forces  in  the  modeled  events,  but  tension  are  introduced  in  both 

Jodels  by  the  inaccurate  hydrodynamic  drag  coefficients  that 
use. 

'■  Material  damping  caused  by  elastic  hysteresis  in  the  cable 
mate»1a^  is  found  to  play  a  significant  part  in  smoothing  the 
compWati on  of  cable  tension,  enen  though  it  has  little  overt 
relation  to  the  overall  event.  Small  amounts  of  material  damp-  , 
ing  ajre- sufficient  to  prevent  undersirable  oscillations  between  I 
model  nodes;  larger  values  enable  SEADYN  to  accurately  model  the 
forced  oscillation  of  the  anchor  as  the  length  of  the  suspension 
dable  was  varied  through  resonance! 

“‘The  laboratory  experiments  were  restricted  to  serially 
connected  cables  and  special  buoys  )r  anchors  suspended  in  a 
vertical  plane  in  still  water  in  order  to  make  the  data  useful 
for  comparing  simple  computer  models  as  well  as  stances  by  com¬ 
paring  equilibrium  and  steady-state  conditions  to  results  ob¬ 
tained  from  elementary  theory.  The  data  from  all  four  experi¬ 
ments  are  available  for  evaluating  other  models. 
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SECTION  I 


INTRODUCTION 


1.1  BACKGROUND. 

Large  rigid  structures  are  supported  by  foundations  resting  deep  in  the 
earth.  This  obvious  convenience  is  manifest  in  most  land  based  civil  engineering. 
When  the  environment  restricts  access  to  suitable  foundation  supports,  other 
design  concepts  must  be  used.  Suspension  bridges  are  one  method  to  simplify 
construction  when  foundation  supports  are  difficult  to  access:  concentrate  the 
rigid  structure  at  a  faw  points.  In  the  ocean,  not  only  the  remoteness  of  the 
foundation,  but  also  the  large  overturning  moments  produced  by  the  drag  from  ocean 
current  militate  against  rigid  structures. 

Cable  structures  in  the  ocean  provide  a  compliant  response  to  hydrodynamic 
forces,  varying  their  geometry  to  provide  more  "leverage"  against  the  flow.  In 
addition,  the  flexibility  of  cables  facilitates  undersea  "construction"  by  compact 
storage  on  shipboard  winches  and  deployment  through  sheaves.  Since  cable  must  be 
loaded  in  tension,  gravity  must  be  counteracted  by  a  support  at  the  top  of  the 
structure.  This  supporting  force  is  nearly  always  provided  by  hydrostatic  buoy¬ 
ancy.  Anchors  represent  the  vestigial  remains  of  the  sophisticated  foundations 
used  to  support  rigid  structures. 

When  the  purpose  of  the  structure  includes  the  need  for  accurate,  consistent 
relative  position  of  parts,  the  cost  of  rigid  structures  in  the  ocean  may  be  ac¬ 
cepted,  as  in  the  offshore  oil  industry,  but  in  deep  water,  cable  systems  are  vir¬ 
tually  always  used. 

Sophisticated  computer  models  have  been  developed  to  analyze  the  stress  and 
deflection  of  rigid  structures.  These  models  can  use  the  nearly  linear  response 
and  small  deflection  inherent  in  most  rigid  structural  materials  to  simplify  the 
solution  of  large  algebraic  matrices.  Computer  models  of  ocean  cable  systems, 
however,  must  contend  with  the  large  deflections  that  these  structures  exhibit,  as 
well  as  the  nonlinear  response  of  many  cable  materials.  Only  in  recent  years  have 
the  mathematical  and  procedural  difficulties  posed  by  the  solution  of  the  response 
of  even  simple  cable  geometries  to  the  dynamic  ocean  environment  been  successfully 
addressed  (Reference  1). 

The  Civil  Engineering  Laboratory  (CEL),  under  the  sponsorship  of  the  Naval 
Facilities  Engineering  Command,  has  engaged  in  a  major  program  to  develop  tech¬ 
niques  for  the  static  and  dynamic  analysis  of  oceanic  cable  structures.  This 
program  is  divided,  on  the  one  hand,  into  a  series  of  small-  and  large-scale  ex¬ 
periments  to  measure  the  response  of  cable  structures  to  their  environment,  and  on 


1-1 


other  hand,  to  develop  and  evaluate  computer  programs  which  attempt  to  predict 
th.it  response. 

The  series  of  experiments  was  devised  to  provide  data  by  which  computer 
models  of  cable  systems  can  be  evaluated.  The  data  from  these  experiments  provide 
not  only  a  standard  for  measuring  the  accuracy  of  model  predictions,  but  also  a 
severe  test  of  the  "robustness"  of  the  computer  algorithms  used  by  the  model. 

Separate  reports  have  already  been  issued  describing  each  experiment  in  de¬ 
tail;  other  reports  compare  computer  model  results  with  the  experimental  data. 
These  reports  are  listed  in  the  References.  This  report  brings  the  results  of  all 
these  comparisons  together  in  summary  form. 

The  series  of  experiments  was  devised  to  test  various  aspects  of  the  models. 
The  experiments  ranged  in  size  from  a  6-foot  laboratory  scale  to  a  2,500-foot 
prototype  mooring  in  the  ocean.  The  small-size  experiments  should  not  be  consid¬ 
ered  as  scale  models  of  large  cable  systems,  but  as  full-size  events  of  their  own. 
This  is  because  all  the  parameters  of  a  cable  system  cannot  be  scaled  simultane¬ 
ously.  The  parameters  for  each  size  experiment  can  be  specified  independently  in 
the  computer  model.  It  is  valid,  therefore,  to  require  that  a  computer  model 
small -size  events  as  well  as  large  ocean  cable  systems. 

1.2  CABLE  GEOMETRIES. 

The  first  experiment  in  the  series  used  thin  elastic  nylon  or  rubber  strands 
about  6  feet  long  (Reference  2).  The  second  experiment  used  a  60-foot  length  of 
the  rubber  material  (Reference  3).  A  third  experiment  (Reference  4)  used  both  the 
nylon  string  as  well  as  the  rubber  material  in  the  6-  and  60-foot  experiments. 

The  length  of  the  cable  was  a  variable  in  this  experiment;  the  length  was  varied 

between  about  6  and  60  feet  by  a  small  winch.  The  fourth  experiment  in  the  series 

was  a  full  scale  mooring  set  in  2500  feet  of  water  (Reference  5). 

The  experiments  included  several  cable  geometries.  The  sketches  on  Fig¬ 
ure  1-1  show  the  cases  that  were  studied.  The  geometries  are  simple,  consisting 
of  serial  cable  systems  with  no  Y  connections,  in  order  to  simplify  computer 
modeling.  The  presumption  is  that  more  computer  programs  will  be  able  to  model 
simple  geometries,  and  that  the  causes  of  node!  aberrations  from  the  data  will  be 
more  easily  identified  if  the  geometry  is  simple.  Nevertheless,  these  cases 
exhibit  many  of  the  dynamic  processes  that  affect  actual  cable  structures  in 
the  ocean. 

1.2.1  THE  PROPPED  PENDULUM.  Figure  l-l(a)  is  a  sketch  of  the  dropoed  pendulum. 
One  end  of  a  cable  is  fixed  at  or  near  the  water  surface.  A  weight  is  attached  to 

the  other  end  and  supported  at  approximately  the  same  elevation  as  the  fixed  end, 

but  displaced  from  it,  so  that  the  cable  hangs  in  a  catenary  arc.  When  the  cable 
is  released,  it  drops,  jerks  the  cable  taut,  then  swings  under  the  fixed  suspen¬ 
sion.  At  the  end,  the  cable  hangs  vertically. 
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Geometry  of  Experiments 


Some  of  the  events  a  computer  model  should  replicate  for  this  geometry  are 
the  static,  elastic  catenaries  that  describe  the  inital  and  final  states  of  the 
cable,  the  "gooseneck"  that  forms  in  the  cable  just  after  the  weight  is  dropped, 
the  occurence  of  slack  cable  at  the  vertex  of  the  initial  catenary  shortly  after 
the  weight  is  dropped,  the  impact  of  the  weight  as  the  cable  comes  taut,  and  the 
damping  of  the  swing  towards  the  final  vertical  shape. 

Although  the  final  configuration  of  this  case  suggests  that  it  represents  a 
measurement  system  suspended  from  a  drifting  surface  float,  these  systems  are  not 
deployed  with  such  brusque  abandon.  It  is  instead  an  approximation  of  the  tech¬ 
nique  used  to  deploy  the  mooring  sketched  in  Figure  l-l(b). 

1.2.2  ANCHOR  LAST  MOORING  DEPLOYMENT.  In  this  technique,  the  buoy  that  supports 
a  mooring  leg  is  towed  in  the  water  behind  a  ship  by  means  of  the  mooring  wire. 

As  the  ship  steams  over  the  desired  site,  the  anchor  is  dropped  into  the  water. 

It  falls  to  the  bottom  dragging  the  buoy  behind  it.  If  the  mooring  wire  is 
shorter  than  the  water  depth  and  the  anchor  is  heavy  enough,  the  buoy  is  sub¬ 
merged.  Figure  l-l(b)  is  a  sketch  of  a  mooring  deployed  by  the  author-last 
technique. 

Computer  models  of  this  case  should  depict  the  initial  towing  "catenary",  the 
horizontal  deceleration  after  the  anchor  is  dropped,  the  "gooseneck",  the  anchor 
speed  and  trajectory,  the  buoy  speed  and  trajectory  (especially  the  growth  of  ten¬ 
sion  that  ultimately  submerges  the  buoy  of  a  subsurface  mooring),  and  the  response 
of  the  mooring  to  the  anchor  impact  as  well  as  the  final  mooring  configuration 
under  the  influence  of  the  prevailing  ocean  currents. 

1.2.3  SINGLE  LEG  MOORING  RELAXATION.  If  a  system  is  mathematically  linear,  then 
its  response  to  an  abrupt  disturbance  is  an  important  measure  of  its  dynamic 
charateristics.  Some  mooring  models  art  based  on  linearized  mathematics,  in  which 
dynamic  effects  are  treated  as  small,  linear  disturbances  about  a  series  of  equi- 
lijrium  states  that  are  updated  non-linearly  every  so  often.  The  mooring  relaxa¬ 
tion  experiments  sketched  in  Figure  l-l(c)  provide  a  step-function  input  to  a 
single  leg  mooring. 

In  the  sketch,  one  end  of  a  cable  is  shown  anchored  to  the  bottom  of  the  test 
tank  or  ocean  with  the  other  end  supported  by  a  subsurface  buoy.  The  buoy  is  dis¬ 
placed  from  its  equilibrium  position,  held  until  equilibrium  is  re-established, 
then  abruptly  released.  The  buoy  floats  back  to  the  original  equilibrium. 

1.2.4  SIMPLE  PAYOUT.  Winches  are  commonly  used  to  deploy  and  recover  cables  at 
sea.  Most  cable  models  to  date  require  fixed  cable  lengths,  but  a  few  can  solve 
the  variable  length  cable  problem.  Figure  l-l(d)  shows  a  weight  being  lowered  or 
raised  on  a  cable. 

1.2.5  SIMPLE  OSCILLATOR.  Figure  l-l(e)  shows  another  comnon  event  when  cables 
are  used  at  sea;  a  weight  is  suspended  from  a  floating  object  that  heaves  in  the 
seaway. 
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1.2.6  PAYOUT  WITH  OSCILLATOR.  This  case  combines  (d)  and  (e)  in  a  more  realistic 
scenariol  using  a  winch  on  a  heaving  ship.  Figure  l-l(f)  shows  the  concept.  The 
experiment  is  deceptively  simple.  Resonant  interactions  between  the  oscillator 
frequency  and  the  varying  cable  length  may  be  masked  in  the  experiment  by  fluid 
and  material  damping,  but  present  an  imposing  obstacle  to  a  computer  algorithm. 

1.2.7  OBLIQUE  OSCILLATOR  WITH  SUSPENDED  LOAD.  Figure  1-1 (g)  shows  a  cable  with 
one  end  fixed  and  the  other  end  constrained  to  move  in  a  vertical  circle.  A  load 
is  attached  at  some  point  along  the  cable.  The  load  may  be  heavy  or  buoyant. 

This  case  represents  the  problem  of  a  ship  anchored  in  a  seaway  or  a  cable  system 
being  installed  from  a  cable  laying  ship. 

1.2.8  TWO  LEG  MOORING  RELAXATION,  TRANSVERSE  DEFLECTION.  Figure  l-l(h)  shows  a 
variation  of  the  relaxation  experiment  described  in  part  1.2.5  above.  The  buoy 
supports  two  mooring  cables  in  this  case,  and  the  initial  deflection  is  along  the 
bisector  of  the  anchor  baseline. 

1.2.9  TWO  LEG  MOORING  RELAXATION,  PARALLEL  DEFLECTION.  The  last  sketch  shows  the 
two  leg  mooring,  but  in  this  case  the  initial  defection  is  in  the  plane  of  the 
mooring.  The  first  case  represents  the  loading  that  results  in  maximum  deflection 
and  holding  power.  This  case  gives  minimum  deflection  per  unit  load,  but  all  the 
load  is  held  by  one  leg. 

Not  every  configuration  shown  on  Figure  1-1  was  tested  at  every  scale. 

Table  1-1  summarizes  the  runs  for  each  experiment. 

The  experimental  data  may  be  used  to  evaluate  many  computer  models  of  cable 
systems.  The  results  from  two  programs  are  described  in  this  report.  SEADYN  is  a 
general  model  for  networks  of  cables  in  three  dimensions,  SNAPLOAO  models  serially 
connected  cables  hanging  in  a  vertical  plane. 

1.3  COMPUTER  MODELS. 

1.3.1  SEADYN.  SEADYN  (Reference  6)  uses  finite  element  techniques  to  analyze 
three  dimensional  networks  of  cables.  The  finite  element  method  divides  the  con¬ 
tinuous  cable  into  segments  (elements)  whose  properties  are  uniform.  Spatially 
varying  forces  are  assumed  to  be  constant  on  an  element.  The  solution  obtained 
for  each  element  is  required  to  align  with  the  solution  obtained  for  adjacent 
elements. 

SEADYN  is  partitioned  into  several  kinds  of  analyses,  including  static  net¬ 
works  (no  acceleration  of  cable  or  fluid),  time  dependent  dynamics,  resonant 
vibrations  (standing  wave  shape),  and  harmonic  analysis  in  the  frequency  domain. 
The  cable  network  may  include  discrete  bodies  such  as  buoys,  anchors,  and  winches 
that  pay  out  or  reel  in  cable.  External  point  loads,  non-steady  three  dimensional 
current  fields,  surface  waves  or  spectra,  and  wind  forces  are  accepted. 

Although  the  finite  element  technique  originated  in  the  analysis  of  linear  or 
nearly  linear  systems,  the  SEADYN  algorithm  accepts  the  nonlinearities  that  result 
from  large  deflection,  large  strain,  fluid  drag,  and  hyper-elastic  materials  with 
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TABLE  1-1.  SUMMARY  OF  EXPERIMENTAL  CASES 
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hysteretic  damping.  The  internal  stress  distributions  produced  by  bending  and 
torsion  are  not  modeled. 

1.3.2  SNAPL.OAD.  SNAPLOAD  (Reference  7)  is  a  lumped  parameter  model  of  serial 
cable  systems  that  hang  in  a  vertical  plane.  This  implies  that  wind,  current,  and 
wave  forces  are  also  aligned  in  that  plane.  The  lumped  parameter  method  replaces 
the  continuous  cable  and  its  partial  differential  equations  with  a  series  of 
"lumps"  where  all  forces  act,  connected  by  massless  springs.  The  algebraic  equa¬ 
tions  of  motion  for  the  "lumps"  are  then  solved.  If  the  parameters  of  the  lumps 
are  properly  chosen,  the  response  of  the  lumps  converges  to  the  response  of 
corresponding  points  of  Lhe  cable. 

Inasmuch  as  serially  connected  two  dimensional  cable  systems  is  the  largest 
class  of  applications,  SNAPLOAD  has  broad  applicability.  Like  SEADYN,  SNAPLOAD 
provides  for  momentarily  slack  cables,  as  well  as  variable  length  cables.  Envir- 
oiimental  forces  from  surface  waves  and  depth-dependent  ocean  currents  are  ac¬ 
cepted.  For  the  class  of  problems  that  it  solves,  SNAPLOAD  is  usually  more 
economical  than  SEADYN. 
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EXPERIMENTS 


2.1  THE  SIX  FOOT  EXPERIMENTS. 

2.1.1  CONFIGURATIONS.  Five  simple  types  of  experiments  were  made. 

2. 1.1.1  Single  Point  Mooring  Relaxation.  In  single  point  mooring  relaxation,  one 
end  of  a  cable  is  fixed  to  the  tank  floor  and  a  buoy  is  attached  to  the  free  end. 

The  buoy  is  displaced  a  specified  distance  horizontally  and  then  released,  allowing 
it  to  return  to  the  original,  or  equilibrium  position,  as  illustrated  in  h'g- 

ure  2-1.  The  variable  for  this  kind  of  experiment  is  the  initial  position  (verti¬ 
cal  and  horizontal  displacement)  of  the  buoy  relative  to  the  fixed  end  of  the  cable 
prior  to  release;  this  detennines  the  initial  configuration  and  tensions  (degree  of 
slackness)  in  the  system.  Twenty-six  of  the  experiments  were  of  this  type. 

2. 1.1. 2  Simulated  Anchor-Last  Deployment.  In  simulated  anchor-last  deployment, 
the  fixed  end  of  tne  cable  is  positioned  just  below  the  surface  of  the  water  and 
an  anchor  is  attached  to  the  free  end.  The  anchor  is  raised  to  a  height  equal  to 
the  fixed  end  and  displaced  a  specific  distance  horizontally,  then  allowed  to 
fall,  and  return  to  the  vertical  equilibrium  position  as  shown  in  Figure  2-2. 

This  experiment  simulates  the  free-fall  anchor-last  deployment  procedure  frequently 
used  for  single  point  moors.  The  test  variable  is  the  horizontal  separation  of 
the  fixed  end  of  the  cable  and  the  anchor.  Sixteen  simulated  anchor-last  experi¬ 
ments  were  made. 

2. 1.1. 3  Bi-moor  Relaxation.  In  a  bi-moor  relaxation,  one  end  of  each  of  two 
cables  is  attached  to  a  common  point  on  a  buoy;  the  other  ends  are  fixed  to  separate 
positions  on  the  tank  floor  as  shown  in  Figure  2-3.  As  with  the  previous  types  of 
experiments,  the  test  variable  is  the  spatial  position  of  the  buoy  prior  to  re¬ 
lease.  For  this  test,  the  buoy  was  displaced  from  equilibrium  either  in  the  ver¬ 
tical  plane  passing  through  the  anchors  or  perpendicular  to  this  plane,  and 
allowed  to  relax  to  the  original  or  equilibrium  position.  The  latter  provided  the 
only  three-dimensional  experiments  in  this  series.  Nine  in-plane  (two- 
dimensional)  and  seven  out-of-plane  experiments  were  made.  The  horizontal  separa¬ 
tion  of  the  fixed  cable  ends  was  5.5  feet. 

2. 1.1. 4  Single  Cable  Suspended  Load.  In  the  single  cable  suspended  load  tests, 
an  anchor  is  suspended  from  a  single  cable  fixed  to  a  load  measuring  device  and 
vertically  oscillating  point  as  shown  in  Figure  2-4.  The  frequency  and  amplitude 
of  the  oscillating  end  of  the  cable  were  varied  to  provide  a  laboratory  simulation 
of  wave-induced  ship  or  platform  motion.  Seven  suspended  load  experiments  were 
made. 
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Figure  2-1.  Diagram  of  Simulated  Anchor-Last  Deployment  Experiment 
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Figure  2- 


.  Diagram  of  Single  Point  Mooring  Relaxation  Experiment 
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Figure  2-3.  Diagram  of  Bi-Moor  Relaxation  Expi 
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2. 1.1. 5  Dual  Cable  Suspended  Load.  In  dual  cable  suspended  load,  one  end  of  each 
of  two  cables  is  attached  to  opposite  sides  of  an  anchor  or  buoy  and  the  remaining 
ends  are  attached  one  to  an  oscillator,  the  other  to  a  load  measuring  device  fixed 
on  the  tank  floor  as  shown  in  Figure  2-5.  The  vertical  and  horizontal  separation 
of  the  cable  ends  was  4.0  and  3.5  feet.  The  oscillating  end  is  moved  in  a  circular 
motion  at  a  predetermined  amplitude  and  frequency.  Two  of  these  dual  cable  tests 
were  performed. 

2.1.2  MATERIALS.  Three  types  of  cable  were  used  in  these  experiments,  a  nylon 
line  and  two  silicon  rubber  cables  used  to  absorb  shock  in  sonobuoys^.  The  nylon 
line  was  inelastic^  whereas  the  rubber  cables  were  considerably  more  elastic.  One 
of  the  sonobuoy  cables  was  solid  rubber.  The  other  included  two  fine  copper  wires 
wound  as  shown  in  Figure  2-6.  The  latter  cable  is  identical  to  that  used  in  the 
60-foot  experiments  and  was  tested  in  the  present  experiments  to  provide  an  assess¬ 
ment  of  scale  effects.  The  physical  properties  of  all  of  the  cables  are  given  in 
Table  2-1  and  the  elastic  behavior  is  plotted  in  Figures  2-7,  2-8,  and  2.9.  All 
the  cables  exhibit  significant  hysteresis  in  their  elastic  properties.  The  data 
shown  were  measured  after  the  cables  had  been  stretched  to  15  percent  of  their 
original  length  five  times  each;  care  was  taken  during  the  experiments  to  insure 
that  the  cables  were  not  stretched  beyond  this  initial  prestretch.  The  black 
rubber  cables  were  marked  at  approximately  6-inch  intervals  with  contrasting  bands 
of  white  silicon  ink;  the  ink  did  not  change  any  of  the  cable  properties. 

The  buoys  and  anchors  used  in  these  experiments  were  fabricated  from  2-inch 
diameter  hollow  plastic  spheres^.  The  surfaces  of  the  spheres  were  smoothed  and 
painted  for  the  experiments.  The  hollow  spheres  were  used  directly  as  buoys;  the 
anchors  were  made  by  lining  the  interior  of  the  spheres  concentrically  with  lead 
shot  to  the  desired  total  weight.  Thus  the  diameters  and  surface  conditions  (and 
drag  coefficients)  of  all  the  buoys  and  anchors  were  the  same.  The  physical  pro¬ 
perties  of  the  buoys  and  anchors  are  presented  in  Table  2-1. 

Complete  assemblies  of  buoys  or  anchors  and  cables  were  made  up  before  the 
experiments  were  begun.  All  of  the  cables  were  mounted  on  the  tension  transducers 
by  a  threaded  ferrule;  a  pin  in  each  ferrule  provided  an  index  for  installation  on 
the  force  transducers  to  give  a  common  zero  point  for  the  length  measurements.  The 
distance  from  the  center  of  the  pin  to  the  center  of  the  buoys  or  anchors  was  1.3 
i nches . 

2.1.3  EXPERIMENTAL  ARRANGEMENT.  The  experiments  were  conducted  in  the  Controlled 
Atmosphere  Launching  Tank  of  the  California  Institute  of  Technology  (CIT)  Hydrody¬ 
namics  Laboratory  during  October  1976.  The  tank  is  a  horizontal  cylinder  13  feet 
in  diameter  and  30  feet  long  (Figure  2-10).  Camera  ports  are  precisely  located 


^  Marsh  Industries,  Mt.  Clemens,  Michigan  48403 

^  Inasmuch  as  all  real  materials  deflect  under  load,  no  material  is  literally 
"inelastic".  The  word  is  used  in  this  report  in  a  relative  sense.  Materials 
that  deflect  slightly  under  a  force  are  inelastic  when  compared  to  materials 
that  deflect  greatly  under  the  same  force. 

3  Emerson  and  Cuming,  Canton,  Massachusetts  02021 
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TABLE  2-1.  MATERIAL  PROPERTIES  FOR  LABORATORY  EXPERIMENTS 


ITEM 

TEST 

SERIES 

DIAMETER 

(in) 

WEIGHT 

IN  AIR 
(lb) 

WEIGHT  IN 
FRESH  WATER 
(lb) 

Smal 1  Buoy 

6  ft 

2.0 

0.025 

-0.121 

Large  Buoy 

60  ft 

2.9 

0.047 

-0.406 

Light  Anchor 

6  ft,  60  ft 

2.0 

0.246 

0.108 

Heavy  Anchor 

6  ft,  60  ft 

2.0 

0.398 

0.252 

Suspended  Load 

Variable 

2.0 

1.16 

1.003 

Sil icon 

0.163 

(Ib/ft) 

8.979x10'^ 

(Ib/ft) 

1.075x10*^ 

Silicon,  wired 
Sample  1 

60  ft 

0.175 

12.8  X  10"^ 

2.00x10'^ 

Sample  2 

60  ft 

0.168 

12.8  X  10"^ 

2.40x10'^ 

Sample  3 

6  ft 

0.163 

11.4  X  10"^ 

2.38x10"^ 

Sil  icon 

Variable 

0.159 

9.2  X  10"^ 

5.82x10"^ 

Nylon 

6  ft 

0.10 

3.73  X  10"^6 

1.057x10"^ 

Nylon 

Variable 

0.059 

.3 

0.792  X  10 

^.73x10 
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Figure  2-7.  Elastic  Characteristics  of  Nylon  Cord 
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Figure  2-9.  Elastic  Characteristics  of  Silicon  Rubber  Cable  with  Conductors 


Figure  2-10,  Six-Foot  Experiment  Test  Facility 


along  one  side  and  at  both  ends.  The  water  in  the  tank  is  treated  and  filtered 
for  optical  clarity.  The  water  in  the  tank  was  8  feet  deep  for  these  experiments. 

Tensions  at  fixed  ends  of  the  cables  were  measured  by  Kistler-Morse  Model  113 
triaxial  strain  gage  deflection  sensors.  Deflection  along  any  axis  produces  a 
voltage  propoKional  to  the  force  component  along  the  axis.  The  sensors  were 
calibrated  before  and  after  the  experiments  and  were  found  to  be  linear  within 
0.05  percent  of  full  scale,  with  no  drift  of  the  instruments  observed.  The  cross - 
axis  interaction  of  the  sensors  was  also  measured  and  found  to  be  less  than 
1  percent. 

Position  time  histories  of  points  on  the  cables  were  measured  by  three 
pulsed-frame  instrumentation  cameras.  The  cameras  were  modified  for  synchronized 
operation  (on-off  and  frame  rate)  for  these  experiments.  The  fields  of  view  for 
the  cameras  were  arranged  to  overlap  in  the  area  in  the  tank  where  the  experiments 
were  conducted.  The  cameras  were  operated  without  shutters  at  10  frames  per 
second.  Xenon  quartz  flash  lamps  keyed  to  the  frame  rate  provided  an  effective 
shutter  speed  of  200  microseconds  at  precise  intervals. 

A  special  mechanism  was  used  to  release  the  buoy  or  anchor  in  the  relaxation 
and  anchor-last  experiments.  In  operation,  the  sphere  is  held  on  the  mechanism  by 
a  slight  vacuum.  At  the  specified  release  time,  the  vacuum  line  is  switched  to  a 
low-pressure  air  source  to  relieve  the  suction  holding  the  sphere.  The  air  pres¬ 
sure  was  carefully  adjusted  to  avoid  giving  the  sphere  an  initial  impulse.  The 
release  mechanism  was  mounted  on  a  pan-and-tilt  head;  for  each  run,  the  angle  of 
the  sphere  was  adjusted  to  miminize  the  initial  rotation  when  the  sphere  was  re¬ 
leased  (Figure  2-11). 

A  variable  frequency  and  amplitude  oscillator  was  used  to  simulate  ship  mo¬ 
tion.  The  frequency  was  controlled  by  changing  the  speed  of  the  driving  motor  and 
was  set  between  0.16  and  1.7  Hz.  The  amplitude  was  adjusted  to  either  1,  3,  or  6 
inches  by  varying  the  radius  of  the  crank  pin  in  a  Scotch  Yoke  mechanism,  as  shown 
in  Figure  2-4.  The  cable  attachment  was  free  to  swivel,  so  that  only  tension  was 
induced  at  the  oscillator  end  of  the  cable. 

Figure  2-12  is  a  photograph  of  a  buoy  relaxation  experiment  rigged  in  the 
empty  tank.  Figure  2-13  shows  a  dual-cable  suspended  load  set-up  using  the  oscil¬ 
lator  with  the  Scotch  Yoke  mechanism  removed. 

The  experiment  was  controlled  through  a  CIT  Graduate  Aeronautical  Labora¬ 
tories'  minicomputer  system.  The  computer  synchronized  the  cameras,  flash  lamps, 
triaxial  force  transducers  and  buoy/anchor  release  mechanism  and  provided  a  common 
time  base  for  all  the  data  and  events  during  an  experiment.  The  computer  also 
controlled  data  acquisition  and  displayed  raw  or  processed  data  for  "quick  looks" 
during  the  experiments. 

Tension  data  were  recorded  for  all  of  the  experiments.  Photographic  position 
data  were  recorded  for  the  relaxation  and  anchor-last  tests. 
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2.1.4  DATA  ACQUISITION.  Photographic  data  were  recorded  at  10  frames  per  second. 
Each  channel  of  the  tri axial  force  gages  was  sampled  and  digitized  every  10  micro¬ 
seconds;  ten  samples  were  averaged  and  recorded  every  0.1  second.  The  data  were 
averaged  to  reduce  noise,  while  the  sampling  rate  was  chosen  to  avoid  aliasing  in 
the  oscillator  experiments. 

2.1.5  DATA  REDUCTION.  The  photographic  data  from  these  experiments  were  reduced 
by  the  Computer  Services  Branch  of  the  Naval  Air  Test  Center,  Patuxent,  MD.  The 
cable  shape  on  each  frame  was  digitized  visually  using  a  film  reader  which  produced 
punched  cards  of  the  results.  The  three-dimensional  nodal  positions  were  computed 
from  the  cards  using  photogrammetry. 

The  nodal  positions  were  merged  with  the  tension  data  at  CEL.  The  zero-time 
reference  was  determined  visually  from  the  photographic  data  by  observing  the 
initial  separation  of  the  sphere  from  the  release  mechanism.  Velocity  and  accelera¬ 
tion  at  a  node  were  estimated  by  differentiating  a  parabolic  function  fit  to  three 
consecutive  position  measurements. 

2.1.6  RESULTS.  The  data  collected  in  this  experiment  are  suitable  for  comparison 
to  numerical  models.  This  is  true  for  both  the  tension  and  photographic  data. 

The  accuracy  of  the  data  in  the  tension  and  nodal  displacement  histories  are 
dependent  on  the  errors  in  the  measured  variables  and  the  time.  In  both  cases,  the 
largest  source  of  error  is  in  the  time  estimate  which  is  limited  by  the  camera 
sampling  period  (0.1  second).  The  actual  time-zero  could  occur  anywhere  within  the 
sampling  period,  with  a  corresponding  shift  in  time  for  the  remaining  data.  This 
undetermined  shift  could  cause  significant  differences  between  the  experimental  re¬ 
sults  and  numerical  predictions  for  the  first  few  points.  As  time  increased  the 
relative  error  in  the  time  estimate  would  decrease. 

The  accuracy  of  the  tension  measurements  is  considered  very  good.  This  judge¬ 
ment  is  supported  by  comparing  the  results  of  calibrations  performed  before  and 
after  the  experiments. 

The  errors  in  the  calculated  nodal  displacements  are  estimated  to  be  approxi¬ 
mately  1/2  inch  in  the  X,  Y,  and  Z  directions.  The  velocities  and  accelerations 
should  be  used  as  estimates  only,  since  they  were  numerically  derived. 

2.2  THE  60  FOOT  EXPERIMENTS. 

2.2.1  CONFIGURATIONS.  Two  simple  types  of  experiments  were  made: 

a.  "Anchor-Last,"  in  which  the  fixed  end  of  the  cable  is  positioned 
near  the  surface  of  the  water  and  an  anchor  is  attached  to  the  free 
end.  These  are  described  in  Section  2. 1.1. 2  and  shown  in  Fig¬ 
ure  2-1.  This  experiment  simulates  the  anchor-last  deployment  pro¬ 
cedure  frequently  used  for  single  point  moors. 

"Buoy  Relaxations,"  in  which  one  end  of  a  cable  is  fixed  to  the  tank 
floor  and  a  buoy  is  attached  to  the  free  end.  These  cases  repeat, 
on  a  larger  scale,  the  6- foot  experiment  described  in  Section  2. 1.1.1 
and  sketched  on  Figure  2-2. 
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2.2.2  MATERIALS.  Both  types  of  experiments  were  made  with  simple  cable  systems, 
with  only  a  single  cable  and  a  single  buoy  or  anchor.  The  experiments  were  also 
two-dimensional,  with  little  or  no  out-of-plane  motion.  The  intent  in  these  ex¬ 
periments  was  to  minimize  the  complexity  of  the  situation  to  be  simulated  by  the 
computer  models,  while  retaining  the  essential  features  of  the  dynamic  behavior  of 
cable  systems. 


2.2.3  EXPERIMENTAL  ARRANGEMENT.  The  experiments  were  conducted  in  fresh  water  in 
the  Hydroballistics  Tank  at  the  Naval  Surface  Weapons  Center  (NSWC),  White  Oak, 
Silver  Spring,  Maryland.  The  tank  is  100  feet  long  x  35  feet  wide  x  75  feet  deep 
and  has  windows  on  the  walls  and  ends  of  the  tank  for  observation  or  photography. 
The  water  depth  in  the  tank  for  these  tests  was  65  feet.  The  water  is  filtered 
and  treated  for  optical  clarity. 

The  models  consisted  of  0.17-inch  diameter  silicon  rubber  shock  mitigation 
cord,  a  2.9-inch  diameter  hollow  plastic  buoy,  and  2.0-inch  diameter  hollow  plastic 
spheres  concentrically  weighted  with  lead  shot  for  anchors.  Properties  of  these 
models  are  presented  in  Table  2-1  and  Figure  2-9.  The  silicon  rubber  cord  had 
two  contrahel ical ly  wound  fine  wire  conductors  which  provided  power  for  small 
light  bulbs  spaced  at  approximately  6-foot  intervals  along  the  cord. 

For  the  6- foot  experiments,  the  cable  was  illuminated  by  high  speed  Xenon 
flash  lights  and  photographed  one  frame  per  flash.  The  60-foot  experiments  were 
conducted  in  darkness.  A  single  photographic  plate  was  exposed  for  each  run,  the 
camera  shutter  opened  at  the  start  of  the  run  and  not  closed  until  the  run  was 
complete.  Once  a  second  during  the  run  the  small  lights  wired  to  the  cable  were 
flashed.  The  images  of  these  lights  define  the  shape  of  the  cable  at  those 
instants. 

^or  each  run  a  buoy  (or  anchor)  was  attached  to  the  free  end  of  the  cable. 

The  fixed  end  was  attached  to  a  Kistler-Morse  Model  113  triaxial  deflection  sensor 
as  described  for  the  6- foot  experiment.  The  sensor  was  calibrated  before  and 
after  the  experiments  and  was  found  to  be  within  0.05  percent  of  the  full-scale 
used.  The  cross-axis  interaction  of  the  sensor  was  measured  and  found  to  be 
within  1 ,0  percent. 

Photographic  coverage  of  the  tests  was  obtained  by  running  the  tests  at  one 
end  of  the  100-foot  tank  and  placing  cameras  at  the  opposite  end.  Three  cameras 
were  used;  one  provided  a  close-up  view  of  the  initial  release  of  the  buoy  or 
anchor,  one  showed  half  coverage  in  the  central  area,  and  the  other  showed  full 
coverage  of  the  cable  motion.  Camera  and  test  apparatus  positioning  are  shown  in 
Figure  2-14. 

At  the  end  of  the  tank,  opposite  the  cameras  and  directly  behind  the  test 
area,  small  lights  were  placed  in  five  of  the  view  ports.  These  light  positions 
provided  references  on  the  photographic  plates  for  use  in  constructing  grids  for 
photo  data  reduction. 

The  Kistler-Morse  triaxial  sensors  were  wired  directly  to  an  analog-to- 
digital  converter;  x,  y,  and  z  forces  were  displayed  on  a  strip  chart  and  digi¬ 
tally  recorded. 
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SIDE  VIEW 


Figure  2-H.  Camera  Positions  for  60  Foot  Experiments 
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2.2.4  DATA  ACQUISITION.  All  recording  equipment  was  turned  on  and  running  and 
the  pulsing  of  the  lights  commenced  prior  to  opening  the  camera  shutters.  When 
the  camera  shutters  were  opened  the  lights  continued  to  pulse  in  their  initial 
positions  for  a  few  seconds  prior  to  buoy/anchor  release.  The  cable  system  was 
then  released  and  allowed  to  return  to  an  equilibrium  position.  The  camera 
shutters  were  then  closed,  providing  a  time  lapse  photograph  of  the  dynamic  se¬ 
quence  on  a  single  plate. 

2.2.5  DATA  REDUCTION.  All  of  the  data  photographs  were  enlarged  to  8-  x  10-inch 
size.  A  photograph  of  the  five  known  position  lights  was  also  enlarged  to  8-x  10- 
inch  size  and  used  to  create  a  data  reduction  grid.  The  grid  was  then  aligned 
over  the  reference  lights  on  each  data  photograph.  The  position  of  the  flashing 
lights  was  read  from  the  grid  and  punched  on  computer  cards  manually. 

The  triaxial  force  data  were  transferred  from  magnetic  tape  to  punched  com¬ 
puter  cards  mechanically.  The  position  cards  and  triaxial  data  cards  were  then 
entered  into  one  computer  program  and  the  data  reduced  to  position,  velocity, 
acceleration,  and  tension  vs  time  to  be  presented  both  in  graphic  and  tabular 
foim.  The  velocity  and  acceleration  were  calculated  from  the  ineasured  position¬ 
time  data. 

During  the  data  reduction  process  various  problems  became  evident.  The  cable 
rotated  slowly  about  its  axis  while  moving  through  the  water,  at  times  obscuring 
the  pulse  lights  from  the  camera.  The  one-second  pulse  was  slower  than  optimum 
during  high  velocity  displacement  and  faster  than  optimum  during  slow  velocities. 
The  blanks  in  the  photographs,  due  to  obscuring  of  the  lights,  and  the  running 
together  of  the  pulse,  due  to  low  velocities,  were  filled  in  manually  where  the 
separation  was  clear  enough  for  reasonable  estimates.  In  other  portions  the  data 
were  omitted. 

2.2.6  RESULTS.  Fifteen  cases  were  run.  Each  run  began  and  ended  with  the  moor¬ 
ing  line  in  a  static  condition.  Since  the  static  shape  and  tension  distribution 
of  a  flexible,  elastic  line  can  be  expressed  by  relatively  simple  formulae,  the 
quality  of  the  dynamic  data  gathered  during  a  run  can  be  inferred  by  comparing  the 
initial  and  terminal  measurements  of  cable  shape  and  tension  to  the  static  calcu¬ 
lation.  Several  runs  were  examined.  Of  these,  three  were  selected  for  analysis  by 
SEAOYN  and  SNAPLOAO. 

2. 2. 6.1  Run  6:  Anchor-Last.  The  positions  of  seven  points  along  the  cable  were 
recorded  prior  to  releasing  the  anchor  and  compared  to  the  static  elastic  catenary 
formulae.  The  mean  radius  from  a  measured  point  to  the  corresponding  point  along 
the  catenary  was  0.5  feet  with  a  standard  deviation  of  0.4  feet.  The  mean  radial 
error  is  about  1  percent  of  the  length  of  the  cable. 

The  cable  force  on  the  load  cell  was  calculated  to  be  0.068  pound  at  81.8  de¬ 
grees  below  horizontal.  The  output  of  the  load  cell  was  0.058  pound  at  an  angle 
of  76.7  degrees.  The  standard  deviation  of  the  load  cell  calibration  is  0.007 
pound. 
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.I;!!  remainder  of  the  discrepancy  must  be  accounted  in  the  standard  deviations  of 
the  ...easured  cable  weight  and  elasticity  used  in  the  catenary  formulae. 

By  90  seconds  after  the  anchor  was  released,  the  mooring  had  virtually  stabi¬ 
lized:  drag  for  :es  were  negligible.  The  elastic  catenary  equations  were  solved 
using  the  measu  ed  weight  and  horizontal  displacement  of  the  anchor  as  boundary 
conditions  for  .his  case.  The  mean  distance  between  the  measured  and  calculated 
node  positions  at  t  =  90  seconds  was  0.7  feet  with  a  standard  deviation  of  0.5 
feet.  This  corroborates  the  1  percent  error  noted  for  t  =  0. 

The  immersed  weight  of  the  cable  and  anchor  in  Run  6  is  about  0.240  pound. 

The  load  cell  output  was  0.221  pound  at  the  end  of  the  run,  100  seconds  after  the 
anchor  was  dropped.  The  load  cell  recorded  the  cable  angle  as  about  5  degrees 
from  vertical.  This  reflects  an  exponential  approach  to  vertical  equilibrium. 

2. 2. 6. 2  Run  11:  Buoy  Relaxation.  At  the  start  of  Run  11,  the  mean  radius  be¬ 
tween  the  calculated  and  measured  locations  of  seven  points  along  the  cable  was 
0.8  feet,  with  a  standard  deviation  of  0.4  feet.  The  initial  tension  was  calcu¬ 
lated  to  be  0.027  pound  at  12.9  degrees  below  horizontal.  The  load  cell  output 
was  only  0.0013  pound  at  37.7  degrees  below  horizontal.  The  tension  after  100 
seconds  was  measured  as  0.271  pound  at  9.2  degrees  from  vertical.  The  tension  due 
to  weight  and  buoyancy  was  0.308  pound. 

The  measured  state  of  this  mooring  at  t  =  90  seconds  was  compared  with  the 
state  calculated  for  a  static,  elastic  catenary,  using  the  buoyancy  and  displace¬ 
ment  of  the  buoy  as  boundary  conditions.  The  mean  radius  from  a  measured  mode 
position  to  its  theoretical  position  was  3.3  feet,  but  the  standard  deviation  of 
the  position  difference  was  only  0.2  feet.  These  values  are  not  consistent  with 
corresponding  values  .rom  Runs  6  and  15.  Run  11  was  accepted  for  comparison  be¬ 
cause  the  nodes  follow  a  smooth  trajecto'y  through  the  relaxation.  The  question 
will  be  addressed  again  in  Section  4. 

2. 2. 6. 3  Run  15:  Large  Anchor-Last.  The  seven  points  compared  for  this  run  aver¬ 
aged  0.8  feet  difference  between  measured  and  calculated  location  prior  to  anchor 
release.  The  standard  deviation  was  0.3  feet. 

The  comparison  at  the  end  of  Run  15  is  as  good  or  better.  The  mean  distance 
is  0.6  feet  with  a  standard  deviation  for  eight  nodes  of  0.4  feet. 

The  calculated  initial  tension  was  0.068  pound  at  82.0  degrees  from  horizon¬ 
tal.  The  load  cell  read  0.060  pound  at  80.1  degrees  from  vertical.  The  combined 
cable  and  anchor  weight  for  this  case  is  0.384  pound. 

2.3  VARIABLE  LENGTH  EXPE.xIMENTS. 

2.3.1  CONFIGURATIONS.  The  variable  length  experiments  model  the  effects  of 
handling  loads  suspended  from  a  ship  in  a  seaway. 

Simple  payout  involves  lowering  a  load  on  an  elastic  line  deployed  from  a 
fixed  winch.  Included  in  payout  is  the  closely  related  case  of  simple  recovery. 
Simple  oscillation  has  the  load  suspended  on  a  fixed  length  of  line  from  a  winch 
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that  is  heaving  in  simple  harmonic  motion.  Both  the  6-foot  and  60-foot  experiments 
included  a  few  cases  with  simple  oscillation.  The  var'’able  length  experiment 
combined  simple  payout  on  reel  in  with  oscillation. 

2.3.2  MATERIALS.  Two  model  assemblies  were  tested.  One  used  an  elastic  cable  to 
support  a  smal 1  anchor,  and  the  other  used  an  inelastic  cable  to  support  the  same 
anchor. 

2. 3. 2.1  Anchor.  The  model  anchor  was  a  2-inch  sphere  weighing  one  pound  in 
water.  The  outside  of  the  sphere  was  painted  in  a  black  and  white  pattern.  The 
cables  ended  in  a  swivel  which  clipped  into  a  stud  on  the  anchor. 

2.3. 2. 2  Cables.  The  dynamic  response  of  a  compliant  system  was  simulated  through 
the  use  of  a  soft  rubber  cable  attached  to  the  anchor.  The  mechanical  properties 
of  this  cable  are  included  in  Table  2-1.  A  plot  of  the  elastic  characteristics  is 
provided  on  Figure  2-15.  The  response  of  an  inelastic  system  was  simulated  through 
the  use  of  a  braided  nylon  line  attached  to  the  anchor.  This  cable  had  approxi¬ 
mately  four  times  the  axial  stiffness  of  the  rubber  cable.  Mechanical  properties 
of  this  cable  are  listfd  on  Table  2-1.  The  results  of  load-deflection  tests  con¬ 
ducted  on  a  6-foot  specimen  of  nylon  line  are  shown  as  Figure  2-7. 

2.3.3  EXPERIMENTAL  ARRANGEMENT.  The  series  of  tests  was  conducted  in  the  High 
Speed  Water  Tunnel  of  the  Graduate  Aeronautical  Laboratories  at  CIT  in  late  1977. 
The  water  tunnel  pumps  were  not  used  during  the  exercise. 

One  of  the  vertical  sections  of  the  tunnel  thus  formed  a  tank  5  feet  in  dia¬ 
meter  and  65  feet  deep.  A  large  port  in  the  pipe  gave  access  above  the  water  sur¬ 
face.  A  mechanism  for  winching  and  oscillating  the  sample  cables  protruded  through 
the  port,  as  sketched  in  Figure  2-16. 

The  traction  winch  mechanism  outside  the  tank  and  the  right  angle  sheave 
cantilevered  over  the  water  rode  on  vertical  tracks  (Figure  2-17).  A  Scotch  Yoke 
mechanism  (Figure  2-18)  moved  the  entire  assembly  in  simple  harmonic  vertical  mo¬ 
tion  along  the  tracks.  A  strain  gauge  measured  the  cable  tension  around  the  right 
angle  sheave  (Figure  2-19).  Oscillation  amplitudes  of  one  tenth  and  one  inch  were 
available  at  frequencies  continuously  variable  up  to  three  cycles  per  second.  The 
winch  could  reel  in  or  out  at  speeds  up  to  2  feet  per  second.  Its  acceleration  was 
also  controllable  up  to  a  maximum  of  4  feet  per  second-squared. 

2.3.4  DATA  ACQUISITION.  Strip  chart  recorders  gave  real  time  verification  that 
the  apparatus  was  functioning  as  desired.  The  permanent  record  of  the  experimental 
data  was  made  on  analog  magnetic  tape.  Five  parameters  were  measured: 

a.  The  length  of  cable  deployed, 

b.  The  winch  speed, 

c.  The  winch  acceleration, 

d  The  tension  of  the  right-angle  sheave,  and 

e.  The  phase  angle  of  the  oscillator. 
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Figure  2-16.  Variable  Length  Experiment  Arrangement 
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Figure  2-18.  Variable  Length  Oscillator  Drive 


Figure  2-19.  Right  Angle  Sheave  and  Oscillator 
Follower  with  Counter  Weight 


The  cable  was  metered  by  rolling  it  between  a  1-foot  circumference  wheel  and 
the  winch  capstan.  A  tachometer  on  the  capstan  drive  provided  winch  velocity.  Its 
signal  was  differentiated  to  give  the  acceleration  of  the  winch.  A  strain  gauge  on 
the  sheave  support  cantilevered  over  the  water  provided  a  tension  signal.  The 
oscillator  phase  angle  was  measured  by  a  rotary  potentiometer  coupled  to  the  Scotch 
Yoke  crank.  A  voice  log  and  a  ran-in-progress  signal  used  the  remaining  two 
channels  on  the  analog  tape  recorder. 

The  data  acquisition  system  and  cable  oscillator  were  started  about  10  seconds 
before  engaging  the  winch  and  continued  to  record  about  10  seconds  after  the  winch 
was  stopped.  Normally  the  winch  was  started  and  stopped  with  a  constant  accelera¬ 
tion  "ramp  function,"  but  occasionally  the  winch  was  abruptly  braked  to  simulate  an 
emergency  stop.  The  cable  meter  wheel  offset  the  actual  cable  length  by  a  total  of 
7.2  feet.  Of  this  bias,  2.2  feet  account  for  the  horizontal  span  from  the  meter 
wheel  to  the  tension  sheave,  2.0  feet  represent  the  distance  to  the  water,  and 
3.0  feet  represent  the  initial  depth  of  the  "anchor"  sphere. 

In  some  cases,  motion  pictures  of  the  anchor  were  taken  through  a  view  port  in 
the  tank  wall. 

2.3.5  DATA  REDUCTION.  The  analog  tapes  were  played  through  an  analog  to  digital 
converter  and  recoTded  on  digital  tape  at  a  rate  of  1000  samples  per  second.  The 
digical  tape  record  consisted  of  values  corresponding  to  the  sensor  output  vol¬ 
tages.  Calibration  of  the  digital  tape  record  to  engineering  units  is  accomplished 
by  a  second  post  processing  computer  program.  This  program  also  smooths  and  tabu¬ 
lates  the  data  for  time  increments  longer  than  0.001  second.  For  the  data  compar¬ 
isons  in  this  report,  the  smoothing  interval  was  0.01  second. 

2.3.6  RESULTS .  Ninety-seven  cases  were  completed,  of  which  55  used  the  nylon  cord 
and  42  use3~tTre  rubber  line.  Thirty-seven  each  of  payout  and  reel-in  cases  were 
included,  with  22  fixed  length  oscillating  cases,  and  one  purely  static  case  with 
no  excitation  whatever.  Thirteen  of  the  payout  cases,  13  reel-in  cases  had  no 
oscillatory  excitation. 

Time  was  not  recorded  during  the  experiment,  except  on  the  strip  chart  re¬ 
corder.  Slight  drift  in  the  analog  recorded  speed  between  the  test  and  digitizing 
is  negligible  as  is  the  inaccuracy  of  the  time  base  used  to  trigger  the  analog  to 
digital  conversion. 

The  tension  data  lie  within  about  0.01  pound  of  the  local  common  trend.  Occa¬ 
sional  values  may  deviate  by  as  much  as  0.05  pound.  These  points  are  obscure  in 
tabular  listings  or  graphic  plots  unless  an  expanded  scale  is  used.  These  errant 
points  have  negligible  influence  on  the  overall  data.  Computer  programs  used  to 
process  the  data  may  need  extra  statements  to  ignore  these  points. 
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Table  length,  was  measured  to  within  about  Q.l  feet.  The  metering  wheel  rolled 
aloii'j  the  relaxed  length  on  payout,  but  measured  stretched  length,  on  reel-in.  This 
is  an  imperceptible  difference  with  metallic  cables,  but  amounted  to  several  feet 
in  cases  using  the  elastic  rubber  line. 

The  winch  elocity  was  monitored  to  within  2  percent.  The  acceleration, 
electronically  differentiated  from  the  velocity  signal,  is  estimated  to  be  accurate 
to  5  percent.  The  oscillator  phase  angle  is  measured  within  2.5  degrees. 

The  four  cable  excitation  channels  may  also  have  occasional  aberrant  values 
much  larger  than  these  error  estimates. 

2.4  ANCHOR-LAST  DEPLOYMENT  AT  SEA. 

2.4.1  CONFIGURATIONS.  The  dynamic  response  of  a  2, 500- foot  long  single  point  sub- 
surface  scientific  mooring  was  measured  at  sea  during  nine  events  including  the 
anchor-last  deployment  of  the  mooring  itself.  The  term  "scientific"  is  used  to 
Jistinguish  the  mooring  from  a  ship  mooring  which  would  use  much  heavier  compo¬ 
nents.  This  mooring  was  fifth  in  a  series  of  six  moorings  deployed  for  the  Mooring 
Dynamics  Experiment  (MOE)  at  the  Pacific  Missile  Range  Facility,  Barking  Sands, 
Kauai i,  Hawaii  in  October  1976  (Reference  8).  The  sixth  mooring  supported  an  array 
of  current  meters  nearby. 

Figure  2-20  is  a  chart  of  the  MDE  environment.  It  shows  the  location  of  each 
mooring  as  well  as  the  deployment  track  for  the  CEL  mooring.  Figure  2-21  is  a 
sketch  of  the  CEL  mooring. 

The  deployment  proceeded  from  the  top  down.  That  is,  the  main  buoy  at  the  top 
of  the  mooring  was  made  up  to  the  mooring  line  and  towed  afloat  behind  the  USNS  DE 
STEIGUER,  steaming  slowly  along  the  deployment  track.  The  mooring  line  was  de¬ 
ployed,  with  instruments  attached  at  appropriate  intervals,  until  the  entire  moor¬ 
ing  except  the  anchor  was  undertow.  The  anchor,  a  large  clump  of  sandbags,  was 
made  up  to  the  mooring  on  a  chute.  It  was  dumped  overboard  at  the  drop  site  and 
plumneted  to  the  bottom,  eventually  pulling  the  entire  mooring  below  the  surface. 
Later,  the  mooring  was  recovered  by  abandoning  the  anchor  using  acoustically  con¬ 
trolled  release  devices. 

The  mooring  tension  was  recorded  by  instruments  at  four  locations  along  the 
cable.  The  positions  of  "pingers"  attached  near  the  tensiometer  were  tracked  using 
the  hydrophone  array  of  the  missile  range.  At  two  of  the  positions,  water  tempera¬ 
ture  and  hydrostatic  pressure  (depth)  was  also  recorded. 

2.4.2  MATERIALS.  Wire  rope  composed  most  of  the  mooring  line.  A  3/ 16- inch  rope 
(3  X  19  construction)  with  a  polyurethane  jacket  was  used.  Short  lengths  of  3/8- 
inch  galvanized  open-link  chain  were  used  to  insert  the  instrument  packages.  One- 
half-inch  chain  and  33  feet  of  3/4-inch  nylon  line  linked  the  acoustic  release 
assembly  to  the  sand  bag  anchor.  Table  2-2  summarizes  the  material  properties. 
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Figure  2-21.  Schematic  Drawing  of  the  CEL  Mooring 
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2.4.3  EXPERIMENTAL  ARRANGEMENT. 

2.4. 3.1  Mooring  Description.  The  primary  flotation  element  was  a  44- inch  diameter 
sphere  producing  1050-pounds  buoyancy  in  seawater.  A  Xenon  flash  lamp  and  radio 
beacon  were  mounted  on  top  as  recovery  aids.  The  mooring  line  was  3/16-inch  dia¬ 
meter  3  X  19  jacketed  steel  wire  rope  with  an  overall  diameter  of  0.390  inch.  As 
shown  in  Figure  2-21,  short  lengths  of  chain  and  nylon  line  were  used  for  assembly 
convenience  and  to  resist  chafing.  Two  Temperature/Pressure  (T/P)  Recorders,  four 
FVRs  and  five  acoustic  pingers  were  installed  on  the  mooring  as  shown.  Two  38-inch 
aluminum  spheres  provided  847  pounds  emergency  buoyancy  in  the  event  the  primary 
buoyancy  failed.  These  two  relatively  large  spheres  were  used  on  the  CEL  mooring 
in  order  to  simplify  the  mooring  as  much  as  possible  for  numerical  simulation. 

Dual  AMF  acoustic  releases  were  used  to  improve  system  reliability.  The  anchor  on 
the  CEL  mooring  was  made  of  sandbags  lashed  together  and  held  in  a  cargo  net 
(donated  by  PMRF).  A  photograph  of  the  anchor  is  shown  in  Figure  2-22.  The  sand¬ 
bag  anchor  was  used  to  reduce  the  possibility  of  damage  to  the  range  hydrophone 
cables  on  the  sea  floor.  A  1500-foot  length  of  1/4-inch  polypropylene  line  with  a 
small  surface  float  was  attached  to  the  main  buoy.  It  was  used  to  displace  the 
mooring  during  the  relaxation  experiments. 

2.4. 3.2  Instrumentation.  The  mooring  was  instrumented  with  two  T/P  recorders  and 
four  FVRs  developed  by  Charles  Stark  Draper  Laboratories  (CSDL,  Reference  9). 

Five  acoustic  pingers  supplied  and  maintained  by  the  Naval  Torpedo  Station  Detach¬ 
ment,  Hawaii  were  used.  The  position  of  these  pingers  was  tracked  by  the  PMR 
hydrophone  array  and  computer. 

Two  types  of  T/P  recorders  were  used  during  the  MDE,  differing  primarily  in 
the  data  sampling  rate.  Both  types  used  calibrated  thermistors  for  temperature 
measurement  and  strain-gauged  diaphragms  for  pressure  measurements.  The  instru¬ 
ments  are  battery-powered  and  self  recording;  the  data  are  digitized  in  the  in¬ 
strument  and  are  recorded  on  a  data  tape  cassette.  The  low-frequency  T/P  (LFTP) 
instruments  sampled  temperature,  pressure  and  time  every  3.75  seconds,  but  recorded 
only  every  64th  sample  (i.e.,  every  4  minutes).  The  high-frequency  instruments 
(HFTP)  recorded  temperature  every  2.08  seconds  and  pressure  every  C.52  seconds  so 
that  four  pressure  words  are  recorded  with  each  temperature  word.  For  the  CEL 
experiment,  a  burst  lasted  28  minutes,  50.56  seconds,  repeated  every  43  minutes, 
41.44  seconds.  Ranges  and  scaling  of  the  instruments  used  in  the  MDE  are  given  in 
Table  2-3. 

The  Force  Vector  Recorder  (FVR)  is  a  self-contained  motion  and  tension  sensing 
package.  The  battery-powered  instrument  records  digitized  data  on  an  internal  tape 
cassette.  The  FVR  has  six  sensor  channels,  including  three  mutually  orthogonal 
servo  accelerometers,  two  magnetometers  for  azimuth  reference  and  a  strain-guage 
load  cell  for  tension  measurement.  For  the  MDE,  tension  cells  with  a  full-scale 
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Figure  2-22,  The  Sandbag  Anchor 


TABLE  2-2.  MOORING  COMPONENTS  ANO  PHYSICAL  DATA 


"Heavy"  in  Seawater 
Buoyancy  in  Seawater 


2-3.  NOMINAL  RANGES  AND  SCALING  OF  TEMPERATURE-PRESSURE  RECORDERS 
EMPLOYED  IN  MOORING  DYNAMICS  EXPERIMENT 


range  of  3,000  pounds  were  used.  In  addition,  an  internal  clock  provides  a  re¬ 
ference  time.  The  instrument  can  be  used  in  a  burst  or  sample  modej  for  the  MOE 
the  burst  mode  was  used,  with  a  hurst  duration  of  14  minutes,  41.2  seconds  and  a 
burst  repetition  period  of  87  minutes,  22.88  seconds.  The  data  channels  are 
sampled  every  0.52  seconds.  The  instrumentation  electronics  are  housed  in  a  13- 
inch  diameter  pressure  shell  which  weighs  85  pounds  in  air  and  45  pounds  in  water. 
In  use,  the  pressure  shell  is  surrounded  by  a  19- inch  diameter  syntactic  foam 
flotation  shell  which  renders  the  unit  neutrally  buoyant  in  seawater,  with  an  air 
weight  of  130  pounds.  The  mooring  attachments  are  aligned  with  the  centerline  of 
the  unit.  Characteristics  of  the  FVRs  are  summarized  in  Table  2-4. 

The  T/P  and  FVR  instruments  are  switched  on  by  the  removal  of  an  external 
magnet.  All  of  the  instruments  (T/P's  and  FVR's)  were  synchronized  to  a  common 
time  base  by  CSDL. 

Five  acoustic  pingers  were  installed  on  the  mooring  to  give  position  data 
using  the  PMRF  acoustic  tracking  range  facilities.  The  pingers  operated  at  fre¬ 
quencies  between  13  and  24  kHz  and  were  pulsed-coded  for  identification.  The  pulse 
repetition  rate  was  1  pulse  per  second. 

2.4.4  DATA  ACQUISITION.  Two  types  of  experiments  were  made  on  the  CEL  mooring 
during  the  MDE.  The  first  of  these  was  the  anchor-last  deployment  of  the  mooring; 
and  the  other,  termed  a  "relaxation"  experiment,  involved  pulling  the  mooring  to 
one  side,  releasing  it  and  allowing  it  to  return  to  equilibrium.  Only  data  from 
the  deployment  are  presented  in  this  report. 

The  anchor-last  technique  is  a  common  method  for  deploying  single-point 
oceanographic  moorings.  In  this  procedure,  the  following  sequence  of  operations  is 
followed.  First,  the  main  buoy  is  released  with  the  mooring  line  attached.  The 
ship  steams  slowly  away  from  the  buoy,  paying  out  mooring  line,  with  any  instru¬ 
ments  or  distributed  flotation  attached.  The  line  is  paid  out  at  approximately  the 
same  rate  as  the  ship  speed  until  only  the  anchor  remains  on  the  ship.  The  anchor 
is  then  released  and  falls  to  the  sea  floor,  pulling  the  mooring  line  and  buoy  into 
the  equilibrium  configuration.  During  the  anchor  descent,  significant  transient 
forces  act  on  the  mooring  line.  The  present  experiment  was  intended  to  measure 
these  forces  and  the  trajectories  followed  by  the  anchor,  mooring  line  and  buoy 
during  the  deployment. 

The  CEL  mooring  was  set  from  USNS  DE  STEIGUER  (T-AGOR-12)  on  1  November  1976, 
at  22°04'  1.4"  N,  159°53"  33.4"  W.  The  track  followed  by  the  ship  during  the  de¬ 
ployment  is  shown  in  Figure  2-20.  The  points  shown  on  the  ship  track  correspond  to 
the  operations  in  the  deployment  sequence  as  follows: 

ALPHA  6  Main  floation  buoy  released 

BRAVO  6  All  mooring  line  and  instruments  deployed 

(Ship  towing  entire  mooring  with  anchor  on  deck) 
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TABLE  2-4.  FORCE  VECTOR  RECORDER  CHARACTERISTICS 


output  for  optimum  accuracy  and  range  in  any  application. 


CHARLIE  6  Anchor  deployed 

DELTA  6  Anchor  position  on  bottom 


It  will 
the  triangle 
ing  location 


be  noted  that  the  mooring  location  is  approximately  in  the  center  of 
formed  by  three  of  the  range  hydrophones.  The  current  meter  moor- 
is  shown  in  the  upper  left  hand  corner  of  the  figure. 


The  anchor  was  released  at  7:14  a.m.;  it  reached  the  bottom  about  9  minutes 
later.  The  weather  during  the  deployment  was  excellent;  the  wind  and  sea  were 
nearly  calm,  with  a  2-foot  swell. 


2.4.4. 1  Data  Acquisition.  The  data  coverage  for  the  CEL  experiment  is  summarized 

in  Table  This  figure  shows  the  instrument  burst  times  and  the  timing  of  the 

several  events  of  the  experiment.  The  anchor  deployment,  the  relaxation  experi¬ 
ments  and  the  firing  of  the  acoustic  release  were  timed  to  coincide  with  FVR 
bursts.  The  data  recovery  rate  for  the  CEL  experiment  was  exceptionally  high.  The 
FVR's  and  T/P's  provided  a  100  percent  data  recovery  throughout  the  experiment. 

Pinger  ALPHA  was  tracked  during  the  deployment,  but  thereafter  its  signal  was 
lost,  apparently  because  it  was  mounted  so  close  to  the  bottom.  Acoustic  "track" 
was  never  established  on  pinger  BRAVO. 

Pinger  CHARLIE  seemed  to  be  tracked  during  the  experiment,  but  further  pro¬ 
cessing  afterwards  revealed  several  large  gaps  in  the  data.  Pingers  DELTA  and  ECHO 
were  tracked  more  reliably,  with  only  occasional  "holidays"  in  their  records. 

2.4.4. 2  Current  Meter  Data  for  MDE  Experiment  5.  Currents  in  the  MDE  area  were 
measured  by  a  vertical  array  of  WHOl  Vector  Averaging  Current  Meters  (VACM's) 
mounted  on  a  separate  mooring  nearby.  The  current  meter  mooring  location  is  shown 
on  Figure  2-20.  The  data  were  reduced  by  the  Pacific  Marine  Electronics  Laboratory 
(PMEL)  of  NOAA.  Current  data  are  presented  in  Appendix  C  of  Reference  5.  In 
general,  current  speed  was  less  than  0,2  knots. 

2. 4. 4. 3  Temperature/Presure  Recorder  Data.  A  low  frequency  temperature  and  pres¬ 
sure  recorder  was  installed  on  the  CEL  mooring  1179  feet  below  the  buoy  at  a 
nominal  depth  of  1385  feet.  This  device  recorded  a  sample  each  4  minutes  continu¬ 
ously  throughout  the  entire  MDE.  Figures  2-23  and  2-24  show  the  temperature  and 
depth  data  recorded  during  the  CEL  experiment  plotted  for  a  period  from  2  hours 
before  the  first  FVR  burst  until  after  the  mooring  was  recovered.  The  depth  plot. 
Figure  2-24,  clearly  shows  the  four  relaxation  experiments.  The  minimum  depth 
which  would  be  recorded  by  the  instrument  was  1066  feet. 

A  high  frequency  temperature/pressure  recorder  was  installed  on  the  CEL  moor¬ 
ing  33  feet  below  the  main  buoy  at  a  nominal  depth  of  244  feet.  The  burst  repeti¬ 
tion  rate  for  the  HFTP  is  twice  that  for  an  FVR.  The  instrument  records  four 
pressure  readings  in  that  interval.  Figures  2-25  and  2-26  show  the  temperature  and 
depth  measured  by  the  HFTP  during  the  mooring  deployment.  The  record  for  other 
bursts  is  plotted  in  Appendix  D  of  Reference  5. 
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COVERAGE  FOR  MDE  EXPERIMENT  FIVE 


Figure  2-23.  Temperature  Data  Recorded  During  the  CEL  Experiment 


Figure  2-24.  Depth  Data  Recorded  During  the  CEL  Experiment 


Figure  2-25.  Temperature  Measured  by  HFTP  During  the  Mooring  Deployment 


Figure  2-26.  Depth  Measured  by  HFTP  During  the  Mooring  Deployment 


2. 4. 4. 4  Force  Vector  Recorder  Data  for  Deployment.  Tensions  measured  by  the  FVRs 
during  the  deployment  of  the  CEL  mooring  are  plotted  as  Figures  2-27  through  2-30. 
The  times  shown  are  elapsed  time  after  the  start  of  the  first  FVR  burst.  The  fig¬ 
ures  are  in  the  order  in  which  the  instruments  were  installed  on  the  mooring.  The 
first  plot  is  for  FVR  1,  installed  nearest  the  top  of  the  mooring,  while  the  last 
plot  is  for  FVR  4,  nearest  the  anchor.  These  figures  are  reproduced  because  of  the 
exceptional  detail  they  present  of  a  complex  dynamic  cable  event.  The  FVR  loca¬ 
tions  are  shown  on  Figure  2-21. 

2.4.4. 5  Acoustic  Position  Data  for  Deployment.  Five  acoustic  pingers  mounted  on 
the  mooring  were  use^to  measure  the  positions  of  points  on  the  mooring  during  the 
deployment.  The  range  did  not  acquire  a  track  on  pinger  BRAVO;  position  data  for 
the  other  pingers  are  presented  in  Appendix  F  of  Reference  5.  Figure  2-31  shows 
the  trajectories  of  the  four  pingers  starting  about  30  seconds  before  anchor  re¬ 
lease.  Figure  2-32  shows  the  shape  of  the  mooring  at  five  times  during  the  des¬ 
cent,  as  defined  by  pinger  position. 

2.4.5  DATA  REDUCTION.  The  MDE  was  a  major  inter-agency  program  to  measure  the 
behavior  of  several  deep-sea  moorings.  In  addition  to  the  CEL  were  the  partici¬ 
pants  listed  in  Table  2-6.  The  experiment  was  extremely  successful,  with  an  excep¬ 
tionally  high  data  recovery  rate.  Each  of  the  participating  activities  took 
responsibility  to  reduce,  analyze,  and  report  the  data  from  their  individual  ex¬ 
periments  and/or  instruments. 

The  data,  reduced  to  engineering  units,  were  submitted  to  NDBO  for  distribu¬ 
tion  and  storage  on  digital  magnetic  tape.  PMR  reduced  the  acoustic  pinger  track¬ 
ing  data.  CSDL  processed  the  data  cassettes  from  the  FVR's  and  T/P's.  NOAA-PMEL 
reduced  the  current  meter  data. 

2.4.6  RESULTS.  The  data  gathered  during  the  anchor-last  deployment  are  of  remark¬ 
ably  high  quality.  The  current  meter  array  confirmed  that  the  currents  during  the 
entire  deployed  time  of  the  CEL  mooring  were  less  than  0.2  knots,  although  varying 
direction  markedly.  These  slow  currents  deflect  the  taut  subsurface  mooring  only 
slightly. 

The  T/P  recorders  confirm  the  stability  of  the  mooring,  with  deflections 
apparent  only  when  the  mooring  was  displaced  deliberately. 

The  four  FVR  records  show  the  anchor-last  deployment  in  remarkable  detail. 

The  major  events  of  the  deployment  are  clearly  portrayed,  from  before  the  anchor 
was  released  until  the  mooring  has  stabilized.  Each  plot  sustains  the  closest 
scrutiny;  relative  tension  according  to  position  in  the  mooring;  time  of  anchor 
release;  time  of  anchor  impact;  distance  from  ship  excitation  vs  distance  from  buoy 
excitation  under  tow;  intensity  of  anchor  impact,  etc. 

Two  of  the  five  acoustic  pingers  were  tracked  most  of  the  time  the  CEL  mooring 
was  in  the  water.  Another  pinger  could  be  tracked  only  intermittently.  Of  the  two 
pingers  near  the  anchor,  one  could  not  be  tracked  at  all.  The  other  was  tracked 
during  the  anchor  descent,  but  could  not  be  tracked  after  the  impact.  It  was 
probably  just  too  close  to  the  bottom. 
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2-28.  MDE  Experiment  5  FVR  2,  Burst  2:  Deployment  Tension  History 


Figure  2-29.  MDE  Experiment  5  FVR  3,  Burst  2:  Deployment  Tension  History 


Figure  2-30.  MDE  Experiment  5  FVR  4,  Burst  2:  Deployment  Tension  History 


DISPLACEMENT  IN  ctET 


Figure  2-31.  Deployment  Trajectory  for  Fingers 
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TABLE  2-6.  PARTICIPATION  IN  THE  MOORING  DYNAMICS  EXPERIMENT 


[  '  1  PARTICIPANT 

RESPONSIBILITY 

< 

National  Data  Bijy  Office 

funding,  buoy  motion  instrumentation, 
data  archiving  and  dissemination 

Office  of  Naval  Research 

funding,  contract  support 

i 

Woods  Hole  Oceanographic 

Institution 

experiment  coordination  and  direction; 
principal  investigator 

■ 

Charles  Stark  Draper 

Laboratories 

manufacture,  preparation  and  opera 
tion  of  mooring  line  motion  and 
temperature/pressure  instruments;  data 
reduction  from  instruments 

Pacific  Missile  Range  Facility, 
Barking  Sands 

acoustic  and  radar  range  tracking 
facilities,  logistic  support 

Naval  Torpedo  Station  Detach 
[rant,  Hawaii 

acoustic  pinger  preparation  and 
operation 

Naval  Oceanographic  Office 

current  meter  preparation  and  re¬ 
furbishment,  data  reduction 

1 

Civil  Engineering  Laboratory 

funding,  mooring  hardware,  technician 
support 

Naval  Facilities  Engineering 

Coiunand 

funding,  current  maters 

i 

i 

Pacific  Marine  Environmental 
Laboratories,  NOAA 

current  meter  preparation  and  re¬ 
furbishment,  data  reduction 

i 

( 

<  i 

EG&G  Washington  Analytical 

Services  Center,  Inc. 

services,  equipment 

2-52 


SECTION  3 


COMPUTER  MODELS 


3 . 1  INTRODUCriON  TO  SEADYN. 

Thfs  general,  three-^mensional  computer  model  of  the  static  and  dynamic  re¬ 
sponse  of  cable  strticcures  to  the  ocean  environment  uses  the  finite  element  method 
of  analysis  (References  6  and  10).  The  system  of  continuous  cables  is  divided  into 
a  set  of  discrete  elements  joined  at  nodes.  The  equations  of  motion  for  each 
element  are  solved  subject  to  two  conditions.  First,  the  ends  of  all  elements 
attached  to  a  common  nede  must  coincide  at  that  node.  Second,  the  forces  acting  at 
a  node  must  sum  to  zero. 

3.1,1.  OPTIONS  AND  C/^PABILITIES.  SEADYN  is  capable  of  performing  nonlinear  ana¬ 
lyses  of  both  branch  and  series-connected  submerged  cable  structures.  Static,  dy¬ 
namic,  nodal,  and  frequency  domain  analyses  can  be  performed  for  complex  systems 
which  include  buoys,  anchors,  fixed  points,  different  cable  materials,  and  payout 
or  r-eel-in.  An  option  exists  to  analyze  mooring  systems  for  surface  ships.  Load¬ 
ings  may  result  from  point  loads,  non-steady  three-dimensional  current  fields, 
surface  waves  (or  spectra)  and  wind  loadings.  Nonlinearities  arising  from  large 
displacements,  large  strains,  velocity- squared  drag,  and  position-dependent  load¬ 
ings,  hyperelastic  materials  and  constraints  on  the  surface  and  bottom  of  the 
current  field  are  all  treated.  Bending  and  twisting  effects  are  not  modelled. 

The  program  is  a  finite  element  model.  Several  solution  techniques  are 
available.  Three  basic  types  of  analyses  can  be  performed: 

•  Static  shape  with  gravity  loads  only, 

•  Equilibrium  shape  with  current  and  point  loads  as  well  as  gravity, 

•  Transient  dynamic  response  of  a  system  disturbed  from  initial 
equilibrium. 

Additional  options  include: 

•  Calculating  natural  frequencies  and  mode  shapes, 

.  Verifying  a  design  for  the  holding  capacity  of  its  anchors,  buoyancy 
of  its  floats,  and  strengths  of  its  lines,  and, 

•  Calculating  the  drag  coefficient  of  strumming  lines. 
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S-toot;  ."I-',  ''-r'-.or;  and  7-tri  tblo -lerigt^.  Sijch  of  tnis  development  has  been 

d^.'in.ted  at  i-'  .'C  .d’,.:!  .,hii;h  plaiiuos  virtually  ai  1  iiionels  of  cable  dynamics: 

r  on'/iir-,]  i  ng  t'  :  'll;!  iCnr!.  '.h;lo  cliis  probloi'i  ‘nas  not  been  entirely  el  iiai  rated , 
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<pu'i^e  to  '  ci  si  111.  Coefficients  of  drag,  elastic  clamping  and  added  mass  are 
•'xa.'.pies.  fha  cxiipurer  model  can  supply  default  values  for  general  cases. 

■  able  3-1  indie..! i.es  the  default  drag  coefficients  r'or  SEADYN  and  SNAPLOAO.  For 
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The  third  category  of  input  parameters  consists  of  parameters  for  controlling 
r.-'e  cuiputer  nr ograii  itself.  Some  of  those  are  useful  and  necessary.  For  exam- 
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TABLE  3-1.  DEFAULT  DRAG  COEFFICIENT 


Cable  Tangentiau  0.025  0.0  all  Tangential  Diamter  Fron 


:il;  .  of  a  finite  element  system  are  simply  joints;  the  equations  of  motion  are 
solved  for  the  elastic  rods  connecting  the  nodal  joints.  The  classical  analysis 
is  based  on  the  differential  equations  of  an  infinitesimal  element  of  cable. 

Snap  loads  occur  when  a  transient  tension  reduction  produces  a  momentarily 
slack  cable.  Th ^y  are  usually  associated  with  impact  events. 

3.3  ACCURACY  AND  COEFFICIENT.  There  are  three  levels  at  which  a  model  performs. 
First,  just  mentioned,  the  model  iiiust  converge  to  a  solution  according  to  some 
criterion.  Second,  the  solution  obtained  should  reflect  qualitatively  the  effects 
of  the  salient  phenomena  occuring  during  the  event  modeled.  Third,  the  solution 
should  be  numerically  accurate  when  compared  to  the  experimental  measurements. 

It  is  important  that  this  numeric  accuracy  be  obtainable  using  characteristic 
parameters  supplied  either  as  default  values  by  the  program  or  by  the  user  without 
reference  to  experimental  data  for  the  case  at  hand.  (The  purpose  of  computer 
models  is  to  redixe  the  need  for  costly  experiments.)  In  the  case  comparisons 
described  below,  most  parameters  are  selected  by  default  or  on  the  basis  of  an  a^ 
priori  engineering  decision.  In  a  few  cases,  the  parameters  were  adjusted  to 
improve  the  comparison. 

For  example,  the  default  drag  coefficient  of  the  anchor  for  the  MDE  mooring 
was  selected  for  a  spherical  anchor  shape.  This  gave  a  descent  speed  substan¬ 
tially  too  large.  The  drag  coefficient  was  adjusted  to  give  a  terminal  velocity 
equivalent  to  the  measured  descent  of  pinger  ALPHA.  This  procedure  is  justified 
in  order  to  show  that  given  correct  coefficients,  SEAOYN  can  give  good  results. 

In  a  design  environment,  the  engineer  must  be  sensitive  to  the  need  for  pertinent 
design  experiments.  A  "handbook"  value  was  not  an  adequate  representation  for  a 
clump  of  sandbags  in  a  cargo  net.  Although  the  wrong  coefficient  allowed  the 
modeled  anchor  to  fall  about  twice  as  fast  as  its  physica"*  counterpart,  the  net 
effect  was  merely  to  compress  the  duration  of  the  event.  Other  effects  were 
slight. 

Whenever  parameters  have  been  adjusted  to  improve  the  model  "fit"  to  the 
data.  It  is  noted  in  the  text. 

3.4  DIFFICULTIES. 

T}iere  are  three  difficult  areas  that  SEAOYN  shares  with  other  dynamic  cable 
models.  Slack  cable,  inelastic^  cable,  and  varying  cable  length  pose  problems 
that  affect  the  efficiency  and  accuracy  of  computer  models. 

3.4.1  SLACK  CABLE.  The  radius  of  curvature  of  a  suspended  cable  varies  directly 
with  the  tension.  In  most  cases  the  radius  of  curvature  is  so  large  that  bending 
moments,  even  in  large  wire  rope,  can  be  neglected.  Virtually  all  cable  models, 
SEAOYN  included,  assume  that  the  cable  is  flexible,  that  is,  it  does  not  produce 
internal  stresses  to  resist  bending. 


^  See  note  2,  page  2-6. 


A  flexible,  suspended  cable  responds  to  low  tension  by  bending  more  sharply 
in  the  direction  that  causes  tension  to  increase.  In  real  cables,  very  low  ten¬ 
sions  mean  that  bending  moments  are  no  longer  negligible,  and  "zero"  tension  does 
not  produce  a  "kink"  in  the  cable  shape. 

Computer  models  discretize  the  cable  into  small  straight  segments.  When  the 
tension  becomes  so  small  that  the  radius  of  curvature  is  less  than  half  a  segment- 
length,  the  segments  are  unable  to  conform  to  flexible  theory.  The  model  response 
becomes  "jerky",  reflecting  the  difference  between  articulation  and  flexibilty. 

The  problem  is  eased  by  using  shorter  elements. 

SEADYN  provides  two  options  to  approximate  the  behavior  of  slack  cable.  Both 
options  allow  elements  to  be  shorter  than  their  relaxed  length.  One  option  com¬ 
putes  a  compressive  load  for  such  elements,  the  other  ignores  the  stress  due  to 
compressive  strain. 

3.4.2  INELASTIC  CABLE.  The  elastic  rigidity  of  a  cable  has  a  strong  effect 

on  the  efficiency  of  T  computer  model.  Changes  in  tension  propagate  along  a  cable 
at  the  speed  of  sound  in  the  cable,  which  varies  directly  with  its  elastic  rigi¬ 
dity.  In  general,  the  time  step  in  a  model  must  be  less  than  the  time  for  sound 
to  traverse  the  length  of  the  shortest  element.  Inelastic  cables  require  short 
time  steps  because  of  the  high  speed  of  sound  (15,000  ft/sec  in  steel  wire). 

Errors  in  the  position  of  the  ends  of  an  element  produce  errors  in  the  ele¬ 
ment  tension  according  to  the  stress-strain  function  for  the  element.  Inelastic 
elements  produce  large  tension  deviations  for  small  nodal  position  errors.  The 
problem  of  inelastic  cables  is  alleviated  by  using  short  time  steps  and  long  ele¬ 
ments.  The  latter  technique  conflicts  with  short  elements  used  to  avoid  slack 
cable  errors. 

3.4.3  VARIABLE  LENGTH.  Cable  systems  that  include  winches  add  a  third  difficult 
area  for  computer  models.  In  order  to  model  the  increasing  length  of  cable,  an 
element  is  allowed  to  increase  its  length  up  to  some  maximum,  then  a  new  node  is 
inserted,  dividing  the  element  and  creating  a  new  sub-element.  The  new  elanent 
then  grows,  and  so  on. 

The  discretization  imposes  a  single  tension  value  for  each  element,  typically 
the  average  of  the  axial  force  at  each  end.  When  a  new  node  divides  an  element, 
the  resulting  average  tension  is  changed  by  the  weight  of  the  length  transferred 
to  the  new  element. 

Each  node  Insertion  produces  a  small  "jolt"  that  propagates  through  the  sys¬ 
tem.  These  tremors  do  not  materially  affect  the  computed  results  unless  the  cable 
system  itself  is  marginally  stable  (very  near  elastic  resonance  with  inadequate 
damping).  The  node  insertions  may  be  visible  on  a  tension  plot.  The  problem  is 
minimized  by  using  short  elements,  frequent  node  insertions  and  short  time  steps/ 
near  node  insertions,  but  this  adds  to  computation  cost.  The  problem  also  occurs 
when  nodes  are  removed. 

The  three  factors  -  slack  cable,  inelastic  cable,  and  variable  length  - 
should  not  be  viewed  as  errors  or  inadeqtiacies  to  be  corrected,  but  rather  as 
"facts  of  life"  to  be  recognized  and  accommodated.  In  most  cases  the  numerical 
effect  is  negligible.  Failure  to  recognize  them  can  increase  the  modeling  cost 
substantial ly. 
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SECTION  4 


MODEL  COMPARISONS 


4 . 1  SIX-FOOT  EXPERIMENT  COMPARISONS. 

4.1.1  OVERVIEW  OF  COMPARISONS.  A  total  of  17  tests  were  selected  from  the 

67  six  foot  tests  for  comparison  to  the  computer  models.  These  tests  were  se¬ 
lected  as  representative  of  the  full  range  of  experimental  conditions  within  each 
of  the  five  types  of  tests.  In  some  cases,  tests  were  chosen  because  test  results 
exhibited  a  distinctive  behavior  or  phenomenon  not  found  in  other  tests.  Table  4- 
1  shows  the  major  paran.eters  of  the  selected  cases.  Table  4-2  summarizes  the 
physical  parameters  used.  Only  spherical  "anchors"  and  "buoys"  were  used  in  the 
6-foot  experiment  in  order  to  test  the  default  drag  coefficients  used  by  the 
models. 

4.1.2  SINGLE  POINT  MOORING  COMPARISONS.  Five  cases  were  selected  for  comparison 
(Reference  12y!  Test  6  used  a  solid  rubber  cord,  drawn  as  taut  as  its  soft  compo¬ 
sition  allows.  Tests  7,  8,  and  9  used  the  three  cords,  all  starting  from  the  same 
slack  initial  catenary.  Test  19  used  the  rubber  cord  with  embedded  conductors, 
drawn  even  tighter  than  Test  6  relative  to  its  elasticity. 

4. 1.2.1  Shape.  Figure  4-1  shows  the  shape  of  the  cable,  as  measured  and  modeled 
at  several  instants  during  the  relaxation  of  the  mooring.  Similar  plots  for  the 
other  cases  differ  only  in  detail.  The  figure  shows  that  the  normal  drag  coef¬ 
ficients  used  in  both  models  were  too  low,  allowing  the  models  to  lead  the  data 
after  the  initial  second. 

4. 1.2.2  Velocity.  Figure  4-2  shows  the  velocity  history  of  the  buoy,  as  measured 
and  as  modeled  by  SEADYN  and  SNAPLOAD.  The  data  show  that  when  the  buoy  is  re¬ 
leased  it  bobs  up  quickly,  initiating  an  impulse  wave  in  the  cord.  The  resonant 
period  of  the  buoy  and  cable  using  elementary  theory^  is  about  0.4  seconds,  which 
seems  to  be  supported  by  the  velocity  trace  on  Figure  4-2.  During  the  first  1.5 
seconds  of  the  test,  both  models  estimate  excess  buoy  velocity.  For  the  remainder 
of  the  test  the  models  agree  with  the  measured  velocity.  Referring  to  Figure  4-1 
and  comparing  5.1  seconds  with  2.0  seconds  shows  that  the  models  do  not  increase 
their  "lead"  perceptibly  after  about  2  seconds  have  elapsed. 

4. 1.2. 3  Tension.  SNAPLOAO  models  the  gross  build  up  of  the  tension,  but  smooths 
virtually  a11  of  the  detail  (Figure  4-3).  This  is  because  SNAPLOAD  invokes  extra, 
artificial  damping  for  a  short  time  after  the  start  of  a  dynamic  event  in  order  to 
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TABLE  4-2.  BUOY,  ANCHOR,  AND  CABLE  PROPERTIES  FOR  THE 
6-FOOT  EXPERIMENTS  MODELING 


BUOYS 

Diameter 

Weight  in  Air 

Weight  in  Water 

2.0  in. 

0.025  lb 

-0.121  lb 

ANCHORS 

Diameter 

Weight  in  Air 

Weight  in  Water 

0.1  lb  Anchor 

2.0  in. 

0.246  lb 

0.108  lb 

0.25  lb  Anchor 

2.0  in. 

0.398  lb 

0.252  lb 

CABLES 

Diameter 

Weight  in  Air 

Weight  in  Water 

Silicon  w/Wire 

0.163  in. 

11.433  X  10’^  Ib/ft 

2.381  X  10’^  Ib/ft 

Silicon 

0.163  in. 

8.979  X  10"^  Ib/ft 

1.075  X  10"^  Ib/ft 

Nylon 

0.1  in. 

3.736  X  10"^  Ib/ft 

1.057  X  10'^  Ib/ft 

-3 


Figure  4-2.  Buoy  Velocity  Comparison  For  Test  8 


maintain  stable  computation.  It  is  not  apparent  whether  SEADYN  is  modeling 
"real"  ten' ’ on  impulses.  The  aberrations  shown  may  be  artifacts  (inter-nodal 
resonance)  that  material  damping  would  remove.  Figure  4-4  shows  what  can  happen 
to  SEADYN  wnen  material  damping  is  not  used  with  a  "stiffer"  material.  This  is 
t.hs  tension  history  for  Test  9,  which  used  the  nylon  cord.  The  amplitude  of 
"ringing"  is  virtually  as  great  as  the  actual  tension  magnitude.  As  before,  the 
artificial  damping  in  SNAPLOAD  has  smoothed  over  all  the  detail. 

4.1.3  SIMULATED  ANCHOR-LAS F  DEPLOYMENT  COMPARISONS.  Anchor-last  deployment 
simulatTon  is  a  substantiaTTy  more  difficult  model i ng  task.  While  the  cases  end 
like  an  inverted  buoy  relaxation,  the  initial  curvature  is  opposite  the  final 
curvature,  and  the  transition  is  accompnaied  by  quasi- impact  condition?  as  the 
anchor  pulls  the  slack  out  of  the  cable  arc. 

a.  Test  32.  This  test  began  with  a  relatively  taut  initial  shape, 
using  the  "softest"  line  iiwterial .  These  conditions  were  expected 
to  produce  the  best  model  comparisons.  Figures  4-5,  4-6,  and  4-7 
show  the  comtarisons.  Both  models  accurately  portray  the  speed  of 
tile  freely  Tailing  anchor,  but  then  understate  the  cable  drag.  The 
predicted  anchor  motion  is  about  3  inches  per  second  too  fast  during 
the  middle  part  of  the  test.  Both  models  accurately  calculate  the 
acceleration  after  the  first  second.  SEADYN  models  the  tension 
history  fairly  well,  but  StNAPLOAD  misses  the  "impact"  of  the  anchor 
as  the  cable  ccmes  taut. 

Test  35.  This  test  used  the  same  initial  shape,  but  substituted  the 
TneTastic  nylon  line  for  the  soft  rubber.  The  heavy  anchor  was  used 
as  well.  SEADYN  modeled  this  cases  with  the  same  strengths  and 
weaknesses  as  Test  32.  SNAPLOAD,  however,  lost  control  of  the 
tension  (Figure  4-3)  and  this  spoiled  the  run.  Perhaps  this  is 
because  the  resonant  period  for  elastic  motion,  based  on  elementary 
theory,  was  about  half  the  period  for  Test  32. 

Test  39.  In  this  test,  the  components  of  Test  32  were  repeated,  but 
the  rubber  '.vas  hung  in  a  catenary  twice  as  deep  (Figure  4-9).  Both 
models  follow  the  gross  motion  at  roughly  the  correct  speed. 

SNAPLOAD  has  not  handled  the  "gooseneck"  that  fomed  above  the 
anchor  just  after  release  very  well.  SEADYN  did  better  during  that 
2-second  interval.  Velocity  and  tension  histories  are  shown  on 
Figures  4-10  and  4-11. 

d.  Test  43.  The  initial  catenary  was  even  more  slack  for  this  case 

(separation  =  2/3  line  length).  The  heavy  anchor  was  used.  Elemen¬ 
tary  theory  predicts  a  resonant  period  of  0.9  seconds  for  elastic 
motion.  Figure  4-12  shows  the  measured  period  to  be  about  0.75  sec¬ 
onds.  SEADYN  was  closer  to  the  elementary  theory,  but  neither  model 
predicted  the  tension  that  actually  occurred. 
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Figure  4-4.  Tension  Comparison  For  Test  9 
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Figure  4-6.  Anchor  Velocity  Comparison  For  Test  32 
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Figure  4-7.  Tension  Comparison  For  Test  32 
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Figure  4-8.  Tension  Comparison  For  Test  35 
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Figure  4-9.  Displacement  Comparison  For  Test  39 
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Figure  4-10.  Anchor  Velocity  Comparison  For  Test  39 
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Figure  4-11.  Tension  Comparison  For  Test  39 


®*  Test  46.  The  initial  separation  between  the  anchor  and  fixed  end 
was  only  0.58  of  the  length.  This  allowed  the  anchor  to  fall 
freely  over  two-thirds  of  the  total  depth  and  produced  an  abrupt 
tension  spike  when  the  inelastic  nylon  line  snapped  taut.  SNAPLOAO 
was  not  compared  to  this  case.  SEADYN  produced  unrealistic  output. 

To  summarize,  a  total  of  five  simulated  anchor-last  deployments  were  tested. 

As  with  the  single  point  mooring  relaxation  tests,  SEADYN  and  SNAPLOAD  performed 
best  for  initially  taut  anchor-last  deployments  with  soft  cables.  Neither  model 
performed  satisfactorily  for  the  slack  tests. 

The  cause  of  poor  model  performance  in  these  slack  tests  has  not  been  deter¬ 
mined.  These  cases  were  modeled  before  the  material  damping  algorithm  had  been  im¬ 
plemented  in  SEADYN.  That  was  found  to  have  a  pronounced  smoothing  effect  on  the 
tension  model.  SEADYN  computed  a  tension  impulse  that  was  larger  if  the  cable  was 
hung  initially  in  a  slack  bight  (compare  Figures  4-7,  4-11,  and  4-12).  This  suggests 
that  the  default  drag  coefficients  by  SEADYN  are  too  low.  Low  anchor  drag  gives 
excess  terminal  velocity;  slack  initial  shape  gives  more  free-fall  time;  low  normal 
cable  drag  reduces  the  effectiveness  of  straightening  the  "gooseneck"  as  a  smooth 
damper  for  the  anchor  impact. 

It  should  be  noted  that  the  SNAPLOAD  results  for  Test  35  represent  the  best 
results  chosen  from  simulation  runs.  The  initial  SNAPLOAD  model  had  lumped  masses 
in  the  same  locations  as  the  nodes  in  SEADYN;  this  was  the  modeling  technique  fol¬ 
lowed  in  all  the  SNAPLOAD  comparisons.  This  initial  model  did  not  produce  realis¬ 
tic  results  -  the  anchor  quickly  dropped  down  and  in  too  far,  compared  to  the 
experimental  data  and  SEADYN  calculations.  Further  runs  were  made  with  SNAPLOAD 
with  the  final  model  using  two  additional  "half"  lumped  masses  at  the  end  of  the 
cable.  The  changes  did  not  produce  any  significant  improvements.  Also,  for  the 
nylon  cables  in  this  and  other  tests,  SNAPLOAD  could  not  converge  on  the  static 
solution,  which  probably  accounted  for  much  of  the  dynamic  instabilities. 

4.1,4  31-MOOR  RELAXATION  COMPARISONS. 

a.  Test  60.  This  test  was  an  out-of-plane  test,  and  could  not  be  mod¬ 
eled  by  the  two-dimensional  SNAPLOAD.  Figure  4-13,  however,  shows 
the  comparison  between  the  experimental  and  SEADYN  tension  his¬ 
tories.  The  calculated  tensions  from  SEADYN  are  extremely  unstable, 
and  show  no  tendency  to  stabilize  even  after  3.5  seconds.  This 
instability  is  characteristic  when  inelastic  cables  are  modeled 
without  material  damping. 

b.  Test  61.  This  test  is  similar  to  Test  60,  except  that  a  silicon 
rubber  cable  was  used  in  place  of  the  nylon  cable.  SEADYN  tensions 
still  show  an  oscillation  about  the  steady-state  value  even  for  this 
softer  cable  material. 

To  summarize,  only  two  bi-moor  relaxation  comparisons  were  made,  and  both 
were  out-of-plane  tests.  Because  these  tests  were  three-dimensional,  only  SEADYN 
was  compared. 
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Figure  4-13.  Tension  Comparison  For  Test  60 


The  results  were  unexpected.  A  close  fit  between  the  experimental  and 
SEADYN  tensions  was  expected  based  on  the  simplicity  and  symmetry  of  the  test,  and 
the  previous  good  comparisons  for  the  single  point  mooring  relaxation  tests.  This 
was  not  the  case,  however,  as  SEADYN  never  converged  to  a  stable  steady-state 
val ue. 

These  simulations  demonstrate  again  that  material  damping  is  essential  for 
good  model  results,  because  the  program  started  form  its  own  initial  equilibrium 
shapes  and  maintained  a  symmetry  of  tensions  and  nodal  positions  in  both  legs  at 
all  times  in  the  dynamic  analyses.  A  small-time  step  was  used  in  the  dynamic 
analyses  and  was  within  the  proper  bounds,  and  was  not  the  cause  of  any  errors. 
Even  small  amounts  of  material  damping  markedly  reduce  these  tension  pulses. 

4.1.5  SINGLE  CABLE  SUSPENDED  LOAD  COMPARISONS. 

a.  Test  67.  For  this  and  the  remaining  comparisons  only  steady-state 
tensions  are  given.  Experimental  results  were  taken  only  after  a 
long  delay  to  allow  start-up  transients  to  die  out;  likewise,  the 
computer  sinulations  were  taken  only  after  a  delay  of  at  least  three 
periods  of  oscillation.  A  best  fit  comparison  was  found  by  shifting 
the  computer  tensions  in  time  to  best  match  the  experimental  tension 
history.  The  tension  comparison  for  this  test  is  shown  in  Fig¬ 
ure  4-14.  Both  models  calculated  reasonable  tensions  compared  to  the 
actual  values. 

b.  Test  68.  The  tension  compa:*i<;on  for  this  test  is  shown  in  Fig- 

ure  4-15.  SEADYN  tensions  macuh  the  actual  tensions  extremely  well; 
SNAPLOAD  tensions  do  not  show  the  regularity  of  the  experimental 
tensions  in  either  time  or  magnitude. 

c.  Test  94.  Both  models  show  good  agreement  with  the  experimental 
results  as  shown  in  Figure  4-16.  Over  this  time  span  the  results  of 
SEADYN  are  slightly  better  than  those  of  SNAPLOAD. 

The  results  of  these  tests  correlate  strongly  with  the  ratio  of  excitation 
frequency  to  the  resonant  frequency  based  on  elementary  theory: 


Test 

Frequency 

Figure 

No. 

Ratio 

No. 

67 

1.19 

4-14 

68 

1.07 

4-15 

94 

0.37 

4-16 

Frequency  was  the  only  parameter  changed  between  runs  67  and  68.  The  experi¬ 
mental  results  on  Figure  4-14  and  4-15  are  essentially  equal,  save  the  frequency 
difference.  In  each  case,  two  of  the  four  cycles  plotted  are  "clipped"  to  zero 
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Figure  4-14.  Tension  Comparison  For  Test  67 
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Figure  4-15.  Tension  Comparison  For  Test  68 
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Figure  4-16.  Tension  Comparison  For  Test  94 


tension,  and  the  other  two  may  in  fact  be  clipped  was  well.  SEADYN  had  not 
achieved  a  steady  state  on  Figure  4-14;  the  trend  was  to  under-predict  the  ampli¬ 
tude.  SNAPLOAD  had  occasional  anomalous  solutions.  Nearer  resonance,  SEADYN 
models  the  data  well,  as  shown  on  Figure  4-15.  The  cable  was  modeled  with  a 
definite  slack  interval  in  each  cycle.  On  Figure  4-15,  the  proximity  to  resonance 
has  exceeded  the  capacity  of  SNAPLOAO. 

Figure  4-16  is  not  only  far  from  resonance,  but  also  the  0.1-pound  load  re¬ 
placed  the  0.25  pound  sphere.  Even  though  the  excitation  amplitude  was  doubled, 
the  dynamic  response  is  only  half  the  mean  value.  Both  programs  were  able  to 
properly  model  this  case. 

4.1.6  DUAL  CABLE  SUSPENDED  LOAD  COMPARISONS.  Test  133  used  the  small  anchor. 

The  two  rubber  cords  hang  in  an  L  shape.  Since  the  anchor  hangs  from  the  upper 
cord,  the  lower  cord  only  provides  “back  tension".  In  Test  142,  a  buoy  was  used. 
That  inverts  the  equilibrium  catenary. 

Figure  4-17  shows  a  sample  of  the  tension  in  the  lower  link  during  Test  133. 
Intervals  of  constant  cension  indicate  the  anchor  was  swinging  in  an  arc  about  the 
fixed  end.  Changing  tension  shows  the  anchor  moving  radially  about  the  fixed  end. 
Both  SEADYN  and  SNAPLOAD  over- responded  to  this  case.  This  is  because  the  tension 
in  the  lower  end  is  sensitive  to  the  terminal  descent  speed  of  the  anchor.  Both 
models  show  the  lower  link  slack  for  fully  half  the  cycle.  This  is  another  confir¬ 
mation  that  the  default  cable  drag  coefficients  in  the  models  are  too  low. 

Figure  4-18  shows  the  tension  for  Test  142.  The  plot  is  much  smoother  because 
the  buoyancy  force  prevents  the  lower  cable  from  going  slack.  In  addition,  the 
force  required  to  displace  the  buoy  laterally  is  less  than  the  force  required  to 
lift  the  anchor  vertically,  so  that  cyclic  variation  in  tension  is  less.  The 
point-to-point  variation  in  the  test  results  reflects  the  natural  period  of  buoy 
heave  (vertical  "bouncing")  on  the  lower  link.  SEADYN  modeled  this  case  well,  but 
SNAPLOAD  appears  to  be  barely  stable.  This  may  be  because  of  the  high  frequency 
of  the  buoy  heave  (6  Hz). 

4.2  THE  60-F00T  EXPERIMENT. 

Of  the  15  experimental  cases  performed,  three  were  selected  for  comparison 
(Reference  13)  with  SEADYN  and  SNAPLOAD.  Run  6  was  a  simulated  anchor-last  deploy¬ 
ment  using  the  0.1- pound  "light"  anchor.  Run  11  was  a  buoy  relaxation,  and 
Run  15  was  another  anchor-last  simulation,  using  the  heavier  0.25  pound  anchor. 

The  two  models,  SEAOYN  and  SNAPLOAO,  are  compared  to  the  data  from  each  of 
the  three  runs  in  three  ways.  First,  the  cable  shape  in  a  vertical  plane  is 
plotted  for  a  sequence  of  instants  during  the  run.  "Snapshot"  plots  are  the 
easiest  to  comprehend,  because  they  show  the  experiment  as  the  observer  would  see 
it:  a  sequence  of  shapes.  Snapshot  plots  show  whether  the  magnitude  of  dynamic 
forces  are  correctly  modeled.  Each  part  of  the  cable  must  follow  its  trajectory 
with  the  correct  acceleration  and  velocity  in  order  to  maintain  accurate  shapes  at 
subsequent  times. 

The  second  basis  for  evaluating  models  is  to  compare  the  trajectories  of 
selected  points  on  the  cable  with  the  trajectories  measured  for  those  points. 
Trajectory  plots  show  whether  the  direction  of  dynamic  force  during  the  run  is 
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Figure  4-18.  Tension  Comparison  For  Test  142 


correctly.  If  the  dynamic  force  is  misaligned,  a  particle  is  accelerated 
into  in  improper  trajectory. 

The  third  portrayal  compares  the  tension  measured  at  the  fixed  end,  plotted 
as  a  function  of  time  during  the  run,  with  the  corresponding  calculated  values. 

4.2.1  RUN  6.  Figure  4-19  is  a  diagram  comparing  the  experimental  cable  with  the 
physical  models  assumed  for  input  to  SEADYN  and  SNAPLOAD.  Nineteen  nodes  were 
used  for  SEADYN;  seventeen  nodes  is  the  maximum  that  SNAPLOAD  can  accommodate. 

Lamp  4  did  not  operate  during  Run  6.  Figures  4-20  and  4-21  are  images  of  the 
punched  card  data  input  for  two  programs. 

Figures  4-22  and  4-23  show  "snapshots"  of  the  shape  of  the  cable  for  a  se¬ 
quence  of  times  during  Run  6.  The  lamp  locations  for  each  time  are  connected  by  a 
dashed  line  to  help  visualize  the  shape  of  the  cable.  The  actual  shape  of  the 
cable  was  a  smooth  arc,  not  the  articulated  shape  shown.  The  models,  however,  are 
also  represented  as  articulated  shapes,  and  it  should  be  remembered  that  they 
.pproximate  the  smooth  arc  by  a  series  of  straight  links. 

The  numbers  at  the  tips  of  the  curves  show  the  time.  The  experimental  times 
are  circled  for  clarity. 

The  initial  node  locations  computed  by  SEADYN  are  in  close  agreement  with  the 
experimental  data,  but  the  initial  static  shape  computed  by  SNAPLOAD  is  in  sub¬ 
stantial  error.  Two  nodes  near  the  vertex  are  so  far  from  their  correct  locations 
that  the  curvature  of  the  cable  is  reversed.  Both  points  are  at  or  near  the 
region  of  maximum  curvature,  where  the  error  required  to  reverse  the  curvature  is 
greatest. 

Both  plots  show  the  extreme  curvature  of  the  cable  as  it  inverts  the  catenary 
while  being  dragged  down  by  the  anchor.  Both  models  produce  corresponding  shapes. 
This  stage  ends  when  the  cable  comes  taut  about  20  seconds  after  the  release  of 
the  anchor.  Both  SEADYN  and  SNAPLOAO  lag  behind  the  data  during  the  first  20 
seconds,  SEADYN  places  the  anchor  at  t  =  20  seconds,  about  where  it  was  at  t  = 

13.5  secconds  -  a  time  lag  of  about  6.5  seconds.  SNAPLOAD  shows  a  similar  time 
lag,  but  only  about  5  seconds.  Both  models  "catch  up"  to  the  data  by  40  seconds 
after  release.  Both  models  show  erratic  cable  motion  at  the  vertex  of  the  cate¬ 
nary  in  the  interval  from  10  to  16  seconds  after  release,  indicating  slack  cable. 

In  the  second  stage  of  the  run,  in  which  the  anchor  is  suspended  from  the 

taut  cable  and  swings  toward  the  fixed  support,  SEADYN  and  SNAPLOAD  calculate  the 
velocity  of  the  anchor  node  to  be  slightly  faster  than  that  experienced  in  Run  6. 

By  t  =  40  seconds,  the  lag  of  the  first  stage  has  been  cancelled. 

The  model  calculations  lead  the  measured  locations  during  the  second  stage  of 
Run  6,  especially  for  the  portion  of  cable  30  to  50  feet  from  the  fixed  end. 
SNAPLOAO  leads  the  data  somewhat  more  than  SEADYN  on  the  average,  even  though  the 
largest  SNAPLOAD  lead,  11  seconds.  Is  only  a  little  more  than  the  maximum  10- 
second  SEADYN  lead. 

Figures  4-24  and  4-25  show  the  path  taken  by  each  lamp  during  Run  6,  as  well 
as  the  path  calculated  by  SEADYN  and  SNAPLOAD.  The  trajectories  during  stage  two, 
when  the  anchor  is  suspended  from  the  cable,  are  essentially  circular  arcs  about 
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Figure  4-19.  Node  Assignments  For  Run  6 
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Figure  4-20.  Data  Deck  For  SEADYN  Model  of  Run 


a 

z 

«-4 

►-  H 

>•  « 

z  z 

Z  Z  Z  UJ 

K  O 

UJ  UI 

UJ  UJ  UJ  a  UJ 

CJ 

»-i  X 

»-  z 

z  z  z  o  a 

Ui 

^  > 

a  « CD 

a  a  a  z  o 

cn  ^ 

t-4  ^  dC  K 

UJ  Z  UJ 

UJ  UJ  UJ 

Z 

z  z  u  en 

azcz^  z^  a 

a  a  a  a  «h  z 

Z  X  U  O  i-< 

<coaui*^aujoui  t 

1  I  1 

o  a  I  z 

X 

3UJu»-<»-uj»-a.zz 

U.  UJ  z  z  ^  u. 

U.U.U.ZUJUJUUJUJUJUJUJK  a 

UJ 

^  a  «  z 

aoLSdcoaocjA^ 

-:^^uaaaaaaaxoa^ 

a 

o  Z  »-•  0.  »>-  M 

Z303CC0UJZ«4 

«««ZOOOOOOOi-i^ZO 

z 

u  M  H  o  cn  K 

MMAUi^Za^UZ 

zzz«zzzzzzzu.a.uju 

M 

o 

o 

o 

a 

o 

a  a 

> 

o 

o 

o  w 

a 

a 

a  o 

c. 

< 

• 

«  « 

•c 

in 

^  to 

a 

♦ 

M 

to 

CJ 

CM 

«4 

N 

o 

o 

90000^0^9  N 

o 

a  a 

aootoo^o'a 

> 

o 

o 

o  to 

M^CMoroarocM 

a 

c. 

a 

a  o 

rv. 

< 

• 

•  • 

»H  •O'O^roa  ♦ 

'C 

in 

*  M 

*-<  fO 

1  q*  ^  ^  CJ  CJ  <0 

a 

a 

1  1  1  I  i  I  I 

O’ 

w 

ro 

M 

CJ 

CM 

>0 

o 

o 

<o 

> 

o 

o  \n 

O' 

(0 

o 

a 

fv 

a 

a  ^ 

rv 

♦ 

•  « 

'O 

in 

m  «H 

a 

^  a 

o- 

M  < 

CJ 

a 

CA  u 

o 

CM 

X 

• 

in 

• 

o  •r  o 

o  ooooooooo  a 

0^ 

CM 

CJ  o 

oarooa>«^naaav 

O' 

00 

o  o 

tH^^cj'oaaoacia 

a 

N 

o  a 

a  ^  • 

rv 

-0 

•  ♦««o*^09^Na  •  •  • 

'0 

n 

a  o  ^ 

CMCJ^^oiO^'^ 

a 

O’ 

M 

CJ 

ro 

(N 

o 

CM 

«i4 

• 

o  a  o 

o  a 

O’ 

CJ  o 

o  a  ^ 

O' 

a  ^ 

^  a 

i-«  0*  o 

a 

rs  j 

O  N 

a  ^  • 

N 

^  « 

*  • 

•  •  to 

'0 

in  u. 

o  o 

a 

^  UJ 

a 

m  u 

o 

ro 

CJ  X 

CJ 

«*4  u. 

• 

«  o 

o  a  o 

o  a  >  o  o  a  o  a  <0 'O  o  ro 

'O  o 

O  V^O'C'COC^r^Cj*  • 

O' 

m  a 

w  a 

■^aooaw^a*'0  O'oo 

a 

rv 

•  rs 

in  0*  • 

. . 

rv 

'O 

o  • 

♦  ^  a 

»H  •  •  ro  a  ro  ro  CJ 

< 

in 

a  o 

^  • 

1  CJ  r^  »H  CM  CJ  i-» 

a 

'C 

1  1  1  1  1  1  1 

O’ 

rt 

ro 

CJ 

CM 

« 

n  >■•  o  m  o 

O  CJ  o 

o  a 

CJ 

»» cn  * 

•  • 

CJ  o 

o  a  CJ 

O' 

o  <  coo 

o  a 

o  o  o 

a 

IS  -1  o 

N 

♦  N 

a  a  * 

N 

^  »*• 

; 

O  • 

•  •  ro 

'O 

m  Qc 

a  o 

a 

«  O  UJ 

O’ 

n  z  ui 

o 

ro 

CM  U  QC 

CM 

^  z  u. 

o 

o  o 

a  a  a  a  o 

9  a  a 

o’oaocjaooo  «-« 

^  o  • 

•  • 

^  •  CJ  o 

o  9  CJ  n  a  a  ro  ro  a  a  O'  * 

O' 

••  a  1  •  o 

^  o 

o  •  a 

9a9CM0>^aacMv4O9  •• 

a 

a  N  o 

o  c» 

a  •  • 

rv 

Z  'C  «0  UI 

• 

•  •^ro*^99a'Oro**  z 

0) 

2  a  ^ 

o 

a  a  o 

CJCJCJ^v'tv^NCJ  z 

a 

3  «  Z  <9 

^  • 

3 

o 

^  n  3  z 

o 

•J 

ro 

O  CJ  a  M 

O 

CJ 

u  cn 

u 

o  o  o  o  o  o 

oooooooooo 

999999909009990 

o 

CJ  w  ^  a  ^ 

Na^o^cjwra'O 

rsa>9««cjro«a*ONaO'0«^ 

CJ 

ooooooooo*^ 

^  ^  ^ 

vi^^cJCJCJCMCJCJCMCMCJCJroro 

ro 

OOOOOOOOOOOOOOOO 

909900090909000 

o 

eooooooooooooooo 

999000909000090 

o 

4-29 


1 


V£) 

C 

3 

QC 


4- 

C 


0) 

•D 

i 


O 

o 


a. 

«ac 


W 

O 


u 

o 


ra 

4-> 

o 


CM 

I 


a; 

u 

3 

CT> 


J 


Jb.  UU  iU.  OU  4b.  00  bO.  00  5b. 


igure  4-2.'  SNAPLOAD  Snapshots  of  Run  6 
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the  fixed  point.  Both  models  meet  this  simple  criterion.  SEADYN  does  not  allow 
the  silicon  cord  to  stretch  quite  enough.  This  is  attributable  to  the  input 
parameters  more  than  the  computer  model . 

The  most  obvious  discrepancies  between  the  models  and  the  data  on  the  trajec¬ 
tory  plots  for  Run  6  occur  for  nodes  7,  8,  and  9  during  stage  one  (anchor  descent) 
Without  support  at  the  anchor  end,  these  nodes  fall  down  and  away  from  the  anchor 
path.  A  little  later,  after  the  anchor  has  passed  by,  the  cable  gets  tighter  and 
the  nodes  are  jerked  back  towards  the  anchor  path.  The  models  allowed  too  much 
time  for  the  anchor  descent.  This  extra  time  allowed  the  nodes  to  deflect  further 
before  the  anchor  pulled  them  back. 

The  measured  trajectory  of  the  anchor  is  virtually  a  vertical  fall  until  the 
cable  snaps  taut.  But  SEADYN  shows  the  anchor  oscillating  slightly  from  side  to 
side  as  it  falls.*  SNAPLOAD,  on  the  other  hand,  calculates  an  anchor  trajectory 
as  much  as  2.5  feet  inside  the  measured  path.  This  is  sufficient  to  infer  that 
the  resultant  force  acting  on  the  anchor  calculated  by  SNAPLOAD  is  substantially 
misaligned. 

Figures  4-26  and  4-27  show  the  tension  at  the  fixed  node  as  calculated  by 
SEADYN  and  SNAPLOAD  compared  to  the  recorded  values.  The  data  record  begins  with 
the  static  value  (0.06  pound)  which  drops  slightly  when  the  tension  wave  from  the 

anchor  release  arrives.  There  is  a  slight  steady  growth  as  support  for  the  weight 

of  the  lower  half-cable  shifts  to  the  fixed  node.  Then  the  tension  abruptly  in¬ 
creases  as  the  cable  straightens,  comes  taut,  and  stretches  to  stop  the  plummeting 

anchor.  This  abrupt  rise  is  damped  smoothly  into  the  steady  value  that  represents 
the  immersed  weight  of  the  cable  and  anchor. 

The  tension  values  calculated  by  SEADYN  and  plotted  on  Figure  4-26  follow  the 
general  trend  of  the  recorded  tension.  However,  large  oscillations  are  calculated 
during  stage  one.  They  are  attributed  to  inadequate  dissipation  of  elastic  energy 

in  SEADYN  because  material  damping  was  not  modeled.  The  rise  in  tension  is  about 

6.5  seconds  late,  corresponding  to  the  lag  noted  on  the  snapshot  plot.  The  final 

equilibrium  tension  is  higher  than  the  measured  value,  and  still  retains  spurious 


*  These  oscillations  are  an  example  of  the  interaction  of  the  user  with  the 
model.  The  node  spacing  used  to  model  these  cases  is  about  10  percent  of 
the  cable  length;  the  radius  of  the  "gooseneck"  that  forms  behind  the 
anchor  is  somewhat  less.  In  reality,  only  a  small  part  of  the  cable  is 
passing  around  the  gooseneck  at  any  time,  and  the  motion  is  smooth.  In  the 
model  with  large  elements,  however,  the  motion  is  exaggerated  and  irregular 
as  the  links  articulate  around  the  sharp  bend.  The  angular  speed  rises  to 
a  peak  and  falls  back  for  each  element.  The  resulting  centrifugal  force, 
which  acts  on  the  entire  element,  therefore  oscillates  exaggeratedly;  its 
vertical  component  is  restrained  and  results  in  negligible  displacement, 
but  its  horizontal  component  produces  the  anomaly  in  the  anchor  trajectory. 


C 


Bkm  nimo 


4-35 


Figure  4-27.  SNAPLOAD  Tension  History  During  Run  6 


oscillations.  SEAOYN  depends  on  measured  values  for  the  immersed  weight  of  the 
anchor  and  cables.  The  discrepancy  between  the  steady  values  on  Figure  4-26  says 
more  about  the  accu^'acy  of  the  tensiometer  used  in  the  experiment  vis-a-vis  the 
scales  used  to  weigh  the  anchor  and  cable.  Figure  4-27  shows  that  SNAPLOAD 
modeled  the  tension  smoothly  throughout  the  event. 

4.2.2  RUN  11  -  BUOY  RELAXATION.  Figure  4-28  shows  how  SEADYN  and  SNAPLOAO  were 
used  to  approximate  the  physical  geometry  of  the  60-foot  buoy  relaxation  experi¬ 
ment.  The  punch-card  input  decks  are  shov;n  on  Figures  4-29  and  4-30  for  the  two 
programs. 

Run  11  is  a  buoy  deflection  experiment.  The  buoy  is  displaced  to  one  side  of 
the  fixed  anchorage,  then  released.  It  bobs  up  and  drifts  towards  a  position  di¬ 
rectly  above  the  fixed  anchorage.  Figures  4-31  and  4-32  are  the  SEADYN  and  SNAP¬ 
LOAD  snapshot  plots  for  Run  11,  respectively.  This  experiment  is  hydrodynamical ly 
identical  to  stage  two  of  Run  6,  with  the  exception  of  the  inversion  of  the 
hydrostatic  force  on  the  cable.  As  in  Run  6,  both  computer  models  lead  the  data. 
By  the  end  of  the  run,  SEADYN  is  roughly  11  seconds  ahead,  while  SNAPLOAD  is  about 
15  seconds  ahead. 

The  measured  data  present  an  anomaly  at  the  start  of  the  run  that  will  be 
shown  more  clearly  on  the  trajectory  plots:  upon  release,  the  anchor  bobs  up,  as 
expected,  but  jumps  away  from  the  fixed  point.  This  suggests  that  the  release 
mechanism  disturbed  the  buoy. 

Figures  4-33  and  4-34  are  the  trajectory  plots  for  Run  11.  They  show  clearly 
the  anomalous  motion  of  the  buoy  at  the  start  of  the  run.  Otherwise,  the  node 
trajectories  are  a  smooth  progression  from  the  initial  point  to  the  final  point. 
Both  SNAPLOAD  and  SEADYN  calculate  a  final  vertical  shape  in  close  agreement  with 
the  calculated  static  results. 

Figures  4-35  and  4-36  show  the  SEADYN  and  SNAPLOAD  comparisons  to  the  meas¬ 
ured  fixed-end  tension.  Both  models  show  the  tension  rising  abruptly  immediately 
upon  buoy  release.  The  data  show  a  2-second  delay  before  the  tension  rises.  Both 
models  converge  to  a  value  of  tension  consistent  with  the  measured  weight  and 
buoyancy  of  the  components  as  provided  in  the  input  data  for  the  models.  SNAPLOAD 
shows  the  tension  reading  its  final  value  within  15  seconds.  The  plot  of  the 
SEADYN  results  more  nearly  approximates  an  asymptomatic  solution,  although  the 
curve  is  irregular  from  about  15  seconds  through  30  seconds  in  the  run. 

4.2.3  RUN  15  -  HEAVY  ANCHOR  DEPLOYMENT.  Run  15  was  a  repeat  of  Run  6,  with  the 
0.1-pound  "light"  anchor  replaced  by  a  heavier  0.25-pound  sphere  of  the  same  size. 
Figure  4-37  shows  the  node  assignments  for  SEADYN  and  SNAPLOAD.  The  data  input 
cards  are  listed  on  Figures  4-38  and  4-39.  The  models  performed  about  the  same  as 
they  had  for  Run  6,  so  those  comments  will  not  be  repeated  in  detail. 

The  snapshot  plots  for  SEADYN  and  SNAPLOAD  are  presented  as  Figures  4-40  and 
4-41.  Both  models  show  the  anchor  falling  more  slowly  than  measured  in  the  ex¬ 
periment  as  was  shown  for  Run  6.  As  in  Run  6,  the  models  overtake  tne  experiment 
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Figure  4-28.  Node  Assignments  For  Run  11 
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Figure  4-29.  Data  Deck  for  SEADYN  Model  of  Run  11 
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Figure  4-31.  SEADYN  Snapshots  of  Run  11 
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Figure  4>32.  SNAPLOAD  Snapshots  of  Run  11 


igure  4-36.  SNAPLOAD  Tension  History 
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Figure  4-37.  Node  Assignments  For  Run  15 
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Figure  4-40.  SEADYN  Snapshots  of  Run  15 
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during  the  pendulum  stage.  The  anchor  used  in  Run  15  was  heavier  than  the  anchor 
for  Run  6.  Less  time  is  available  for  slack  cable  to  develop  at  the  vertex,  and 
both  models  as  well  as  the  data  reflect  the  smoother  descent.  This  is  also  appar¬ 
ent  in  the  trajectory  plots  (Figures  4-42  and  4-43)  for  SEADYN  AND  SNAPLOAD. 

The  fixed-end  tension  plots  on  Figures  4-44  and  4-45  show  that  both  models 
and  the  data  converged  to  the  immersed  weight  of  the  anchor  and  cord  with  good 
accuracy.  SEADYN  gives  somewhat  erratic  values,  as  in  the  other  cases.  SEADYN 
calculated  the  time  for  the  anchor  to  jerk  the  cord  taut  about  5  seconds  late; 
SNAPLOAD  was  1.5  to  2  seconds  late. 

4.3  THE  MOORING  DYNAMICS  EXPERIMENT  MOORING  SIX. 

Mooring  Six  of  the  Mooring  Dynamics  Experiment  is  sketched  in  Figure  2-21. 

The  deployment  of  this  mooring  by  the  anchor-last  technique  was  modeled  using  the 
SEADYN  computer  program  (Reference  14).  SNAPLOAD  was  not  compared  to  this 
experiment. 

Figure  4-46  is  a  schematic  diagram  of  the  node  model  used  to  approximate  the 
physical  mooring.  Images  of  the  SEADYN  input  cards  are  shown  on  Figure  4-47. 
Appendix  A  reviews  the  development  of  the  parameters  describing  each  node.  Other 
parameters  were  accepted  as  provided  by  default  in  the  SEADYN  model,  with  one  ex¬ 
ception.  The  drag  coefficient  of  the  sandbag  anchor  clump  was  selected  by  match¬ 
ing  the  terminal  velocity  of  the  anchor  to  the  descent  speed  measured  for  pinger 
ALPHA  during  the  deployment. 

4.3.1  SNAPSHOT  GEOMETRY  COMPARISONS.  Figure  4-48  shows  the  shape  of  the  MDE/CEL 
mooring  calculated  at  selected  times  during  the  deployment.  At  time  t  *  0,  the 
anchor  is  on  the  towship  at  the  upper  right  corner  of  the  figure,  and  the  main 
buoy  is  on  the  water  surface  at  the  upper  left  corner  of  the  figure. 

Even  though  the  ship  is  towing  the  mooring  from  left  to  right  at  3.25  knots, 
the  anchor  does  not  coast  an  appreciable  distance  after  release.  Indeed,  it 
swings  back  first,  under  the  recovery  buoys  attached  to  the  acoustic  release. 

Then  it  falls  nearly  vertically  down  the  right  margin  of  the  figures,  until  the 
main  buoy  is  pulled  under  about  6  minutes  (360  seconds)  after  release.  Then  the 
anchor  begins  a  sedate  pendulum- swing  to  the  left  under  the  main  buoy  until  it 
impacts  on  the  bottom  during  the  eighth  minute  after  release. 

The  actual  anchor-last  deployment  differs  from  the  6-  and  60-foot  simulations 
not  only  in  scale.  The  upper  end  of  the  simulated  mooring  was  fixed,  in  order  to 
enable  the  tension  to  be  measured  there.  In  the  at-sea  deployment,  there  is  no 
fixed  point.  The  buoy  is  towed  along  the  water  surface  by  the  descending  anchor, 
using  the  water  as  a  great  "sheave-wheel"  (the  "water-pulley").  This  accounts  for 
the  different  appearance  of  Figure  4-48  from  snapshot  plots  of  the  simulated  de¬ 
ployments,  like  Figure  4-41.  Furthermore,  the  simulated  anchor  weights  did  not 
impact  the  bottom  of  the  test  tank. 

It  is  easy  to  overlook  another  implication  of  the  actual  deployment.  Once 
the  anchor  was  "dropped"  in  the  SEADYN  model,  there  was  no  fixed  point  to  restrain 
the  overall  position  of  the  nodes  in  the  computer  model.  Residual  errors  in  the 
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Figure  4-43.  SNAPLOAD  Trajectories  For  Run  15 
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Figure  4-44,  SEADYN  Tension  History  During  Run  15 
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Figure  4-45.  SNAPLOAD  Tension  History  During  Run  15 
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Figure  4-46.  The  SEADYN  Model  of  MDE  Experiment  Five 


Figure  4-47.  SEADYN  Input  Data  Deck 
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Figure  4-47.  SEADYN  Input  Data  Deck  (Continued) 


TIME  (SECONDS  AFTER  ANCHOR  RELEASE) 


acceleration  or  velocity  could  introduce  increasingly  erroneous  positions  for  the 
system.  Even  though  no  experimental  data  are  plotted  on  Figure  4-4B,  it  connctes 
a  verification  of  the  internal  consistency  of  the  model. 

4.3.2  TRAJECTORY  COMPARISONS.  Figure  4-49  shows  the  paths  taken  by  major  nodes 
during  the  deployment.  The  trajectories  should  be  viewed  individually;  viewed 
collectively,  they  can  be  confusing.  Start  with  the  anchor  at  node  one.  Upon 
release  it  swings  under  the  recovery  buoys,  pulling  them  under  one  after  the 
other,  then  plummets  straight  down  the  right  margin  of  the  plot  until  the  main 
buoy  pulls  under,  when  the  anchor  begins  to  swing  to  the  left.  Finally,  it  im¬ 
pacts  on  the  bottom  at  a  depth  of  2465  feet. 

FVR  4  at  node  three  coasts  a  moment  after  release,  then  is  jerked  under  by 
the  anchor  and  plurrmets  down  with  it. 

FVR  3  at  node  eight,  the  temperature/pressure  (T/P)  sensor  at  node  11,  and 
FVR  2  at  node  13  all  coast  a  little  while,  then  are  pulled  down  in  sweeping  arcs 
following  the  "water  pulley"  principle  (the  drag  on  a  cable  perpendicular  to  it  is 
much  greater  than  the  drag  tangent  to  the  cable.  When  a  cable  is  pulled  sideways 
in  the  water,  it  tends  to  bend  as  if  the  water  were  a  great  pulley  wheel.  The  ca¬ 
ble  tends  to  pull  "around  the  corner"  more  than  it  "cuts  across"  the  corner). 

After  the  anchor  hits  bottom  these  three  nodes  abruptly  stop  falling  and  swing  on 
arcs  about  the  anchor  as  fixed  point. 

FVR  1  at  node  21  and  the  main  buoy  at  node  22  are  towed  along  under  or  on  the 
surface  while  the  rest  of  the  cable  goes  around  the  "water  pulley".  Then  the  an¬ 
chor  pulls  the  main  buoy  under,  and  shortly  thereafter  impacts.  These  nodes  then 
pivot  up  and  to  the  right  about  the  anchor. 

Figure  4-50  is  a  duplicate  of  Figure  4-49  with  locations  of  the  nodes  meas¬ 
ured  during  the  MDE  superimposed.  It  is  easy  to  see  the  triangle  symbols  (pinger 
ECHO)  closely  paralleling  the  trajectory  of  node  21.  The  circle  symbols  (pinger 
ALPHA)  follow  node  three  behind  the  anchor.  The  square  symbols  (oinqer  DELIA) 
loosely  follow  node  eight.  The  hexagons  (oinger  CHARLIE)  clearly  follow  node  13. 

4.3.3.  TENSION  COMPARISONS.  Figures  4-51  through  4-54  show  the  calculated  tension 
superimposed  on  Figures  2-27  through  2-30.  The  resemblance  is  striking,  except 
for  the  large  oscillations  of  tension.  The  large  tension  pulse  produced  by  the 
anchor  impact  is  shown  on  Figur -s  2-27  through  2-30  to  have  decayed  fully  within 
30  seconds;  one  may  infer  that  other  inputs  decay  in  a  like  time.  Therefore,  the 
tension  oscillations  on  those  plots  are  primarily  due  to  continuous  excitation, 
namely  the  action  of  ocean  waves. 

The  SEADYN  model  contained  no  simulation  of  wave  action.  One  would  therefore 
expect  smooth  tension  traces  with  discrete  perturbations  lasting  about  30  seconds 
coincident  with  anchor  release,  anchor  impact,  and  perhaps  the  immersion  of  the 
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Figure  4-49.  Modeled  Trajectories  of  Major  Nodes 


Figure  4-51.  Tension  Comparison  For  FVR  1 


Figure  4-52.  Tension  Comparison  For  FVR  2 


igure  4-53.  Tension  Comparison  Tor  FVR  3 
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main  buoy.  This  is  found  to  be  the  case  for  some  elements  some  of  the  time: 
element  eight  from  about  3.5  to  7  minutes  after  anchor  release,  for  example,  on 
Figure  4-53.  The  large,  undamped  oscillations  of  tension  shown  on  Figures  4-51 
through  4-54  (especially  Figure  4-52)  are  artifacts  of  the  undamped  SEADYN  model. 

These  oscillations  represent  the  resonant  exchange  of  elastic  energy  between 
adjacent  elements.  This  energy  is  dissipated  in  nature  by  material  damping,  i.e., 
hysteresis  in  the  stress-strain  curve  of  the  element  material.  The  SEADYN  algo¬ 
rithm  did  not  include  material  hysteresis  when  this  comparison  was  performed. 

The  orag  coefficient  of  the  anchor  and  buoys  was  calculated  (Appendix  A)  as 
part  of  the  input  model  to  produce  terminal  velocities  equal  to  those  measured  in 

MDE,  so  the  time  from  release  to  impact  is,  not  surprisingly,  nearly  equal  for  the 

model  and  experiment. 

However,  the  beginning  and  duration  of  sub-events  muy  be  usefully  compared, 
after  aligning  the  impact  pulse.  For  example,  both  data  and  model  indicate  that  a 
very  low  tension  occurs  16  to  18  seconds  after  anchor  release.  Figure  2-29  shows 
the  main  buoy  pulled  under  97  seconds  before  impact.  This  is  marked  by  the  change 
from  large,  wave-induced  tension  pulses  to  smaller  double  pulses  after  immersion 
(tension  waves  echoing  off  the  massive  main  buoy).  Buoy  immersion  is  shown  as 
96  seconds  before  impact  on  Figure  4-52  from  the  SEADYN  results. 

The  post- impact  "ringing"  had  a  4-second  period,  as  shown  on  the  FVR  plots  in 
Section  2.  The  corresponding  plots  in  this  section  show  a  period  of  about  4.5 

seconds  for  the  first  few  oscillations  after  impact.  This  "ringing"  dies  out 

after  3-4  cycles  on  the  FVR  plots  and  either  dies  out  or  is  replaced  by  a 
spurious  oscillation  on  the  SEADYN  plots  within  four  cycles  of  impact. 

At  the  start  of  the  anchor  drop  the  vertical  force  exerted  by  the  falling  an¬ 
chor  is  turned  by  the  "water  sheave"  principle  into  a  roughly  horizontal  force  on 
the  main  buoy.  The  tension  is  defined  primarily  by  the  hydrodynamic  drag  tangent 
to  the  cable  and  nodes.  But  as  the  anchor  falls  deeper  and  the  water  sheave  is 
worn  away,  the  cable  tends  more  and  more  towards  the  vertical.  The  main  buoy  is 
pulled  deeper  in  the  water  until  it  comes  awash  and  submerges.  Both  model  and 
data  reflect  this  smooth  increase  in  tension  followed  b/  a  period  of  roughly 
constant  tension  until  the  anchor  impacts  on  the  bottom. 

4.4  VARIABLE  LENGTH  COMPARISONS. 

A  schematic  diagram  of  the  variable  length  experiment  was  sketched  in  Fig¬ 
ure  2-16.  Ninety-seven  runs  were  performed  using  various  combinations  of  ai\ 
oscillation  frequency  and  payout/reel -in  rates.  Two  material  types  were  used,  an 
inelastic  nylon  line  and  a  plastic  rubber  cord. 

The  two  models,  SEADYN  and  SNAPLOAD,  are  compared  to  the  data  by  examining 
the  tension  vs  time  plots  for  each  run.  Tension  was  measured  at  the  top  node 
(payout/oscillation  point).  Visual  scans  of  the  output  were  made  to  check  that 
the  bottom  node  position  appeared  reasonable.  Direct  comparisons  are  not  possible 
since  the  anchor  position  was  not  tracked  during  the  experiment.  Reference  15 
describes  the  comparison  of  SNAPLOAD  and  SEADYN  to  the  variable  length  experiment 
data  in  detail . 
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Appendix  B  lists  the  data  input  decks  for  each  run.  SEADYN  used  a  standard 
node  spacing  of  5  feet  for  all  runs.  The  SNAPLOAD  spacing  was  adjusted  to  accom- 
moaate  modeling  deficiencies  in  the  program.  Very  long  segments  were  used  with 
SfiAPLOAD  in  order  to  avoid  inserting  nodes  during  payout. 

In  the  actual  experiment  the  oscillator  was  about  2  feet  above  the  water  sur¬ 
face.  Neither  program  could  model  this  set-up  exactly  so  the  approximation  was  made 
that  the  entire  cable  was  underwater.  The  difference  in  total  drag  and  weight  of 
the  cable/anchor  system  due  to  this  approximation  is  insignificant. 

Elasticity  tests  gave  two  tension/strain  curves  for  nylon,  one  for  loading  and 
one  for  unloading.  Neither  curve  seemed  appropriate  for  SEADYN  when  the  resonance 
peaks  were  compared  with  both  payout  and  reel -in  runs.  A  stiffness  (EA)  value  of 
133  lb  was  then  computed  for  Run  23  by  assuming  resonance  at  2.5  Hz  and  applying  the 
formula  for  the  natural  period  of  a  mass  on  a  spring.  This  effective  value  gave 
good  results  for  other  runs,  both  payout  and  reel -in,  and  was  within  the  experi- 
imentally  determined  range  of  100  to  200  lb. 

The  nylon  material  damping  was  also  adjusted  to  better  correlate  SEAOYN  with 
the  data.  The  damping  was  increased  to  1.0  lb-sec  from  a  predetermined  value  of 
0.05  lb-sec. 

SNAPLOAD  does  not  model  material  damping.  Mathematical  damping  was  eliminated 
by  setting  the  artifical  damping  factor,  DMAX,  equal  to  0.0001,  and  the  ramp 
limit,  TDE,  to  99999.  The  smooth  results  were  obtained  using  exceptionally  small 
error  limits:  A8EBD  =  RGEBD  *  0.0001.  SNAPLOAD  relies  on  hydrodynamic  damping  to 
limit  resonant  response. 

All  SEAOYN  runs  ran  to  conclusion.  SNAPLOAD  encountered  certain  difficulties 
during  Runs  12,  47  and  48.  Run  12  is  a  reel -in  run  and  it  was  discovered  that 
SNAPLOAD  stops  when  only  two  cable  segments  remain  underwater.  Since  the  run  was 
two-thirds  complete,  it  was  not  rerun. 

During  the  SNAPLOAD  model  of  Run  47,  the  time  step  decreased  to  a  very  small 
value.  The  program  continued  to  execute  until  the  computer  time  allotted  was 
exhausted  but  the  model  time  barely  changed.  SNAPLCAD  does  not  limit  the  number 
of  times  the  time  step  is  decreased,  so  it  is  possible  to  waste  computer  resources 
on  ineffective  computation. 

The  SNAPLOAD  code  does  not  allow  element  strains  greater  than  0.5.  Modeling 
Run  48  aborted  when  the  soft  rubber  exceeded  this  limit. 
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4.4.1  RUN  3.  Run  3  is  a  nylon  payout  without  oscillation.  Figure  4-55  is  a  ten¬ 
sion  vs  time  plot  of  the  SEADYN  model  and  the  data.  Figure  4-56  is  the  correspond¬ 
ing  SitAPLOAD  plot.  The  SNAPLOAD  plot  contains  anamolous  spikes  that  represent 
instantaneous  numerical  instabilities  in  the  program,  perhaps  due  to  too  large  a 
time  step  or  error  tolerance.  SEADYN  produces  a  smooth  curve  for  the  entire  run. 
Previous  experience  using  SEADYN  without  material  damping  resulted  in  tension 
plots  containing  many  undamped  oscillations.  Tiny  steps  are  visible  on  Figure  4-55 
when  SEADYN  inserted  nodes  for  the  payout. 

Both  models  almost  identically  reproduce  the  data.  Both  models  have  a 
slight  tension  offset  of  about  0.05  pound.  This  could  be  caused  by  drag  modeling 
deficiencies  in  the  models,  an  error  in  measuring  payout  rate,  or  simply  a  data 
collection  error  for  the  tension. 

The  first  four  and  final  six  seconds  of  the  run  represent  a  cable  accelera¬ 
tion  and  deceleration.  SEADYN  appears  overdamped  at  the  beginning  and  underdamped 

near  the  end.  SNAPLOAD  exhibits  good  initial  transient  response  but  is  under- 
damped  at  the  end. 

4.4.2  RUN  23.  Run  23  is  a  nylon  payout  with  oscillation.  Figures  4-57  and  4-58 

show  the  corresponding  plots.  The  data  show  an  increase  in  peak  to  peak  tension 

near  the  time  at  which  the  natural  frequency  of  the  cable/anchor  system  approached 

the  oscillating  frequency. 

Both  models  exhibit  this  resonance  near  the  proper  time;  the  SNAPLOAD  irooel 
peak  is  displaced  by  about  1.5  seconds.  SNAPLOAD  is  slightly  underdamped  at  reso¬ 
nance.  The  two  models  produce  almost  identical  transient  responses  at  the  end  of 
the  run.  Both  models  underpredict  this  response,  perhaps  due  to  incorrect  damping 
at  this  point. 

4.4.3  RUN  12.  Run  12  is  a  nylon  reel-in  with  oscillation.  Figures  4-59  and  4-60 
show  the  corresponding  plots.  SNAPLOAD  terminated  at  the  point  where  only  two 
cable  segments  remained  in  the  water.  Up  to  this  point  the  results  of  both 
SEADYN  and  SNAPLOAD  are  in  excellent  agreement  with  the  data.  The  SNAPLOAD  phase 
is  offset  slightly  since  SNAPLOAD  does  not  allow  for  an  initial  phase  shift  in  the 
input  data. 

SEADYN  models  the  transient  portion  at  the  end  of  the  run  well,  picking  up 
the  negative  peak  at  24.5  seconds.  The  SEADYN  response  perhaps  is  slightly  under¬ 
damped  towards  the  end  of  the  run. 

4.4.4  RUN  29.  Run  29  is  a  rubber  payout  with  no  oscillation.  Figures  4-61  and 
4-62  show  the  corresponding  plots.  The  SEADYN  curve  contains  small  discontinui¬ 
ties  that  represent  cable  tension  interpolations  at  node  insertions  (mitosis 
points).  These  "steps",  though  small,  are  much  larger  than  those  noted  on  Fig¬ 
ure  4-55.  An  error  has  subsequently  been  corrected  in  the  SEADYN  code.  Both 
models  follow  the  transient  shape  at  the  beginning  and  end  of  the  run  with  over¬ 
prediction  in  both  models. 
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Tension  Vs  Time  For  Variable  Length  Test  12  (SEAOYN) 
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Figure  4-52.  Tension  Vs  Time  For  Variable  Length  Test  29  (SNAPLGAD) 


3oth  niodels  show  a  tension/time  slope  in  the  constant  velocity  portion  of  the 
run.  The  SNAPLOAD  slope  is  about  twice  that  of  the  SEADYN  slope.  This  slope  is 
caused  oy  the  frictional  drag  on  the  cable;  as  niore  cable  enters  the  water  the 
frictional  drag  increases.  The  SNAPLOAD  frictional  drag  coefficient  is  programmed 
as  a  fixed  0.025  while  SEADYN  uses  a  default  value  of  about  0.018.  Both  values  are 
above  the  normal  range  of  0.002-0.008,  introducing  an  increasing  error  as  cable  is 
paid  out. 

Subsequent  analysis  of  the  nylon  run  showed  the  same  trend,  although  it  was  not 
as  evident  since  the  nylon  cable  diameter  was  less  than  the  rubber  diameter. 

4.4.5  RUN  47.  Run  47  is  a  rubber  payout  with  oscillation.  Figures  4-63  and  4-64 
show  the  corresponding  plots.  Again,  the  effect  of  the  frictional  drag  error  is 
evident.  Discounting  this  error,  both  models  give  a  good  approximation  of  the  data. 
SiiAPLOAO  did  not  complete  the  run  within  the  allotted  time. 

Both  models  appear  to  underdamp  the  initial  transient  and  overdamp  the  final 
transient.  SNAPLOAD  also  overdamps  the  final  half  of  the  constant  velocity  portion. 
No  explanation  was  determined  for  this  difference  in  response.  SEADYN  models  both 
the  high  frequency  excitation  and  the  low,  natural  frequency  after  the  winch  was 
braked.  The  modeling  of  the  high  frequency  tension  amplitude  runs  counter  to  the 
data  when  the  two  frequencies  are  superposed.  The  data  show  an  increase  in  high 
frequency  tension  amplitude  after  the  winch  stops  (from  0.06  lb  before  braking  to 
0.C8  lb  peak  to  peak  after  braking).  SEADYN  is  virtually  exact  before  braking,  bur 
models  the  tension  as  decreasing,  instead  of  increasing,  to  0.03  lb  after  braking. 
Ironically,  this  minor  discrepancy  suggests  slightly  excess  material  damping. 

4.4.6  RUN  48.  Run  48  is  a  rubber  reel-in  with  oscillation.  Figure  4-65  shows  the 
SEADYN  plot  of  this  run.  SNAPLOAD  aborted  Run  48  when  the  st**ain  exceeded  0.5. 
Except  for  the  frictional  drag  error,  the  SEADYN  simulation  is  quite  good.  The 
transients  appear  overdamped  at  the  beginning  and  underdomped  near  the  end,  but  the 
trends  exactly  follow  the  data. 

4.4.7  SUMMARY  OF  VARIABLE  LENGTH  COMPARISONS.  SEADYN  modeled  the  characteristic 
response  of  all  runs  accurately.  Minor  difficulties  pertain  to  damping  in  the 
transient  range,  small  time  step  sizes  required  by  node  mitosis,  and  a  possible 
frictional  drag  coefficient  error.  The  most  striking  results  of  these  runs  is  seen 
when  the  SEADYN  results  are  compared  with  the  plots  of  modeled  tension  for  the 
earlier  experiments.  The  addition  of  material  damping  has  eliminated  the  large 
tension  abberrations  that  marred  the  earlier  modeling,  a  significant  improvement. 
Overdamping  or  underdamping  occurred  frequently  because  the  material  damping 
parameter  was  only  crudely  estimated. 

SNAPLOAD  could  model  some  runs  accurately,  but  the  frictional  drag  coefficient 
error  is  significant.  One  run  stopped  due  to  an  insufficient  number  of  segments 
remaining  in  the  water.  A  rerun  with  additional  segments  presumably  would  work. 
SNAPLOAD  will  not  model  cables  in  which  the  strain  is  greater  than  0.5.  Without 
artificial  or  material  damping,  overdamping  or  underdamping  occur  frequently.  The 
excessive  tangential  cable  drag  coefficient  replaces  some  of  the  lost  damping. 
Tension  spikes  occurred  during  one  run,  possibly  due  to  an  improper  error  tolerance. 
The  most  serious  problem  was  the  failure  of  SNAPLOAD  to  abort  an  ineffective  run. 
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Figure  4-63.  Tension  Vs  Time  For  Variable  Length  Test  47  (SEADYN) 
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Figure  4-64.  Tension  Vs  Time  For  Variable  Length  Test  47  (SNAPLOAD) 


SECTION  5 


OBSERVATIONS  AND  CONCLUSIONS 


5.1  EXPERIMENTS  AND  DATA. 

Four  experiments  have  been  performed  that  measure  the  dynamic  excitation  and 
response  of  cable  systems  in  water.  These  experiments  are  characterized  by  size: 
6,  60,  2,500  feet,  and  variable  length.  The  small  size  experiments  are  not  scale 
models  of  their  large  size  systems;  they  represent  small  systems  in  their  own 
right.  Each  experiment  has  produced  data  suitable  for  evaluating  computer  models 
of  cable  systems.  Simple  shapes  were  used  in  order  that  the  widest  range  of 
computer  programs  would  be  compatible  with  the  geometry.  For  the  same  reason, 
all  the  cables  were  suspended  in  a  vertical  plane. 

The  data  from  each  experiment  series  have  been  shov;n  to  be  rigorous  in  that 
the  computer  models  were  severely  taxed  in  order  to  obtain  solutions.  Some  cases 
are  less  severe  than  others:  buoy  relaxations  and  simple  payout/reel -in  do  not 
demand  as  much  as  the  anchor  drops  or  oscillation  with  payout  or  reel -in  through  a 
resonance. 

Tracking  the  position  of  cable  points  throughout  an  event  was  a  difficult 
measurement.  The  method  used  for  the  6-foot  experiment  was  accurate  to  within 
about  8  percent.  Selected  events  in  the  60-foot  series  were  found  with  accuracies 
estimated  to  be  within  1  or  3  percent  of  the  cable  length.  Misalignment  of  the 
coordinate  grid  during  data  reduction  could  introduce  larger  errors  that  were  hard 
to  detect.  The  Mooring  Dynamics  Experiment  (MDE),  conducted  at  sea,  relied  on 
acoustic  tracking  of  pingers  attached  to  the  mooring.  Positions  could  be  measured 
within  about  6  feet  when  the  tracker  was  "locked"  onto  e  solution.  However,  the 
depth  measurement  could  abruptly  jump  by  several  hundreo  feet  if  the  tracker  be¬ 
came  undecided  among  multiple  solutions. 

Tension  was  the  most  accurate  measurement  in  all  four  experiments.  The 
errors  seemed  greatest  in  the  laboratory  tests  when  the  cables  were  nearly  slack. 
Then  residual  turbulence  in  the  water  and  cable  stiffness  became  perceptible.  The 
direction  of  the  tension  force  was  affected  more  than  the  magnitude,  based  on 
comparing  equilibrium  measurements  from  the  60- foot  experiment  with  the  static, 
elastic  catenary  equations. 

Measurement  of  the  size  and  weight  of  components  was  essentially  routine. 
However,  the  displacement  of  the  model  buoy  and  anchor  spheres  should  not  be 
estimated  from  their  nominal  shape. 


5-1 


The  nylon  and  silicon  rubber  materials  used  for  the  laboratory  experiments 
exhibited  plastic  properties:  hysteresis  in  their  load-deflection  curves.  Creep 
affected  the  rubber  cable  in  the  variable  length  experiments.  Not  only  did  the 
winch  pay-out  cable  in  a  relaxed  state  and  reel  it  in  stretched,  but  the  relaxed 
length  after  the  reel -in  was  longer  than  before  the  deployment. 

Linear  cable  models  require  the  user  to  supply  an  equivalent  elastic  modulus 
for  plastic  cable  materials.  Furthemore,  these  experiments  have  amply  demon¬ 
strated  the  importance  of  internal  material  damping  to  the  stability  of  tension 
calculation  in  numerical  models.  Without  material  damping,  the  calculated  tension 
shows  large  irregular  variations  that  appear  abruptly  in  some  element  and  then 
abruptly  shift  to  another  element.  For  SEADYN,  at  least,  even  unrealistically 
small  values  of  material  damping  are  sufficient  to  eliminate  the  undulations.  But 
larger  (and,  one  hopes,  correct)  values  allow  the  computer  to  portray  the  tension 
accurately  even  when  the  cable  is  oscillated  near  resonance.  This  was  shown  in 
the  variable  length  cases. 

For  most  ocean  cable  systems  the  goal  is  to  avoid  resonances.  The  internal 
damping  of  the  cable  has  negligible  effect  on  the  design.  But  it  is  clear  that  a 
reasonably  accurate  value  is  important  to  the  design  process,  enabling  the  com¬ 
puter  models  to  operate  efficiently  and  portray  their  results  unmasked  by  spurious 
val  ues. 

The  variable  length  experiments  may  be  an  effective  procedure  to  measure 
material  damping  and  equivalent  elastic  modulus.  The  cable  length  for  which 
SEADYN  models  resonance  is  sensitive  to  the  elastic  modulus;  the  amplitude  of  the 
resonant  tension  oscillation  computed  by  SEADYN  is  sensitive  to  the  material 
damping  parameter.  Reference  16  reports  the  current  knowledge  of  the  properties 
of  nylon  line.  Experiments  like  the  variable  length  runs  with  oscillation  could 
be  perfomed  for  larger  sizes  in  a  deep  test  tank  or  lake.  Very  long  lengths 
could  be  tested  at  sea,  but  the  frequency  analysis  would  be  substantially  compli¬ 
cated  by  the  unavoidable  wave  spectrum. 

5.2  THE  SEADYN  MODEL. 

It  should  be  stated  at  the  start  and  perhaps  repeated  several  times  in  the 
discussion  below  that  solutions  obtained  from  the  SEADYN  model  have  consistently 
modeled  the  major  responses  of  cable  systems.  It  is  a  point  easily  stated  and 
easily  forgotten  in  a  discussion  focussed  on  the  errors  of  detail  in  the  model 
comparisons. 

SEADYN  has  been  developed  significantly  during  this  series  of  comparisons. 
Early  efforts  were  devoted  to  the  determination  of  which  solution  algorithm  ap¬ 
plied  to  which  problem.  Then  the  ability  of  the  algorithm  to  converge  to  the 
solution  was  strengthened.  Material  damping  was  added  only  in  time  for  the  last, 
variable  length,  series  of  comparisons;  the  dramatic  improvement  in  the  tension 
model  is  apparent.  Even  small  amounts  of  material  damping  eliminate  the  "ringing" 
that  masked  the  tension  in  earlier  comparisons.  The  variable  length  results 
suggest  that  the  "right"  value  of  damping  enables  SEADYN  to  model  thr  tension 
during  resonant  oscillation  accurately. 

There  are  several  other  areas  of  "unfinished  business"  for  the  SEADYN  model. 
These  do  not  affect  the  underlying  accuracy  of  the  model,  but  improve  its  ability 
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to  model  details  of  dynamic  events,  its  efficiency  in  obtaining  results,  and 
reduce  the  technical  skills  required  to  use  the  program  properly. 

Variable  length  problems  require  the  SEADYN  model  periodically  to  add  or  de¬ 
lete  nodes  by  splitting  or  merging  elements.  This  is  done  to  keep  the  variable 
length  elements  from  becoming  too  long  or  short.  SEADYN  requires  a  very  small 
time  step  during  this  process  of  changing  nodes,  called  mitosis.  The  time  step 
size  control  algorithm  in  SEADYN  does  not  properly  account  for  this  process.  The 
current  approach  is  simply  to  force  the  entire  problem  to  run  with  a  small  time 
step.  At  press  time,  this  problem  had  been  confirmed  by  other  SEADYN  users.  A 
coding  error  was  identified  and  corrected. 

A  third  approach  to  the  mitosis  problem  was  tried  that  may,  with  improvement, 
be  useful.  The  time  step  algorithm  can  efficiently  reduce  the  step  size  for  a  mi¬ 
tosis:  it  halves  the  time  step  iteratively;  but  after  the  mitosis  is  complete, 
the  recovery  to  a  larger  step  is  very  slow.  In  fact,  the  model  takes  as  long  or 
longer  with  this  approach  as  it  does  when  a  small  step  is  specified  for  the  entire 
run.  The  improvement  suggested  is  to  inrrease  the  step  size  faster  after  a  mi¬ 
tosis.  This  requires  that  the  algoritnm  be  able  to  distinguish  mitosis  from  other- 
events  requiring  step  changes.  The  recent  correction  may  obviate  this  approach. 

There  is  a  small  error  in  the  tension  that  is  inherent  to  the  discretization 
of  the  continous  cable  to  finite  elenents.  The  tension  for  an  element  is  calcu¬ 
lated  as  the  average  of  the  forces  applied  at  its  ends.  When  a  node  is  inserted 
or  removed  during  mitosis,  this  average  is  changed  by  the  weight  and  drag  of  the 
length  of  cable  split  off  or  merged.  Step  changes  in  the  tension  are  apparent  in 
the  variable  length  plots,  marking  the  mitosis  events.  In  a  cable  system  of  ocean 
engineering  size,  these  discontinuities  are  negligible.  They  can  be  minimized  by 
using  shorter  elements  if  necessary,  but  there  is  a  cost  penalty  in  execution 
time. 

The  algorithm  that  computes  the  equilibrium  shape  of  a  cable  system  can  find 
the  initial  state  from  very  simple  approximations  given  by  the  user.  The  anchor 
last  deployments  were  defined  with  the  cable  horizontal  and  under  no  tension;  the 
buoy  relaxation  cables  were  input  as  vertical  lines.  This  is  a  convenient  form  to 
define  the  element  lengths. 

The  treatment  of  hydrodynamic  drag  on  the  cable  elements  is  deficient.  The 
default  functions  for  the  drag  perpendicular  to  an  element  consistently  produce 
insufficient  normal  drag.  The  result  is  that  transverse  cable  motions  are  too 
rapid.  The  variable  length  experiment  was  the  only  series  sensitive  to  tangential 
drag.  It  shows  that  the  default  functions  overpredict  these  forces. 

Examination  of  the  SEADYN  code  shows  that  the  default  drag  functions  are  not 
coded  in  the  subroutine  DRAGCO  where  they  would  be  expected;  indeed,  an  undocu¬ 
mented  set  of  functions  occur  there.  The  only  way  to  avoid  the  default  functions 
is  to  write  a  substitute  subroutine  named  DRAGCO  and  re-compile  the  program  with 
it.  This  requires  significant  computer  skills  of  the  user  as  well  as  delays  the 
solution.  The  user  is  required  to  know  the  FORTRAN  computer  language,  understand 
the  concepts  of  subroutine  and  argument,  be  able  to  compile  the  code  and  link  it 
into  the  existing  SEADYN  program.  While  this  procedure  is  a  valid  ultimate  re¬ 
source,  it  should  not  be  required  simply  to  change  a  drag  coefficient. 


These  results  have  shown  that  SEADYN  represents  a  sound  mathematical  analy¬ 
sis.  Significant  progress  has  been  made  in  developing  it  from  a  prototype  in  the 
computer  laboratory  to  a  production  code  available  for  ocean  engineers.  Finally, 
repeating  the  first  observation,  SEADYN  always  displayed  the  characteristic  be¬ 
havior  of  the  cable  system  in  view. 

5.3  THE  SNAPLOAD  MODEL. 

SNAPLOAD  was  used  successfuly  to  model  most  of  the  experimental  cases 
studied.  Each  program  should  be  viewed  within  its  own  context  rather  than  by 
comparison  to  the  results  of  the  other  program,  because  each  program  was  written 
with  different  intentions.  SNAPLOAD  is  an  attempt  to  model  serially-connected 
cables  suspended  in  a  vertical  plane  with  reasonable  accuracy  and  economy.  Inas¬ 
much  as  the  majority  of  ocean  cable  systems  fall  in  this  class,  the  potential 
tradeoff  of  numerical  sophistication  against  computer  economy  and  user  simplicity 
is  attractive.  It  is  “fair",  however,  to  compare  the  actual  cost  of  using  the  two 
programs  in  these  baseline  studies.  The  SEADYN  program  was  steadily  updated  while 
these  cases  were  modeled,  but  very  little  program  maintenance  was  done  on 
SNAPLOAD. 

Although  SNAPLOAD  is  bcsed  on  a  different  mathematical  analysis  than  SEADYN, 
the  program  demonstrated  the  same  difficulty  in  modelling  tension  when  inadequate 
damping  was  present.  The  model  has  a  mathematical  damping  algorithm,  but  no 
formula  is  supplied  to  relate  the  mathematical  damping  parameter  to  cable  proper¬ 
ties.  The  algorithm  was  intended  to  stabilize  the  program  during  the  early 
seconds  of  a  run;  it  decays  along  a  "ramp  function"  whose  duration  is  specified  by 
the  user. 

SNAPLOAD  requires  a  good  initial  estimate  of  the  shape  of  the  cable  system. 
Rough  estimates  require  substantial  computer  time  for  SNAPLOAD  to  refine.  The 
crudest  estimates  (entering  the  nodes  in  a  horizontal  or  vertical  row  for  conven¬ 
ience)  may  not  converge  at  all.  This  is  not  a  major  problem  because  many  programs 
are  available  to  calculate  the  equilibrium  shape  of  serial  cables  in  a  vertical 
plane. 

One  of  the  premises  for  SNAPLOAO  is  that  it  is  easy  to  use.  The  node  number¬ 
ing  scheme  counts  by  four,  and  is  divided  among  "major  nodes",  where  cable  proper¬ 
ties  may  change  and  discrete  loads  inserted,  and  ordinary  nodes  that  simply  divide 
a  cable  into  parts.  A  revision  is  in  process  to  simplify  the  nomenclature  for 
users. 

The  hydrodynamic  drag  coefficients  are  not  specified  by  the  SNAPLOAD  user. 

The  normal  drag  on  cables  was  consistently  understated  in  these  runs,  allowing 
both  anchor-last  deployment  simulations  and  buoy  relaxation  cases  to  translate  a 
little  too  rapidly  transverse  to  the  cable  axis.  The  tangential  drag  coefficient 
was  significant  only  in  the  variable  length  cases.  The  value  coded  in  the  program 
overpredicts  the  tangential  drag  in  all  these  cases. 

SNAPLOAD  will  only  run  until  170  solution  times  have  been  tabulated,  because 
the  solution  tables  are  accumulated  in  the  computer  core  memory.  Tnis  can  be  in¬ 
creased  by  recompiling  the  program  with  larger  DIMENSION  arguments  for  the  arrays 
used.  A  better  approach  involves  recoding  the  program  to  write  the  arrays  to  a 
binary  file. 
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Cable  payout  is  restricted  to  cases  in  which  the  payout  link  can  accommodate 
the  entire  payout  length.  In  the  variable  length  cases,  SNAPLOAD  computed  an 
extremely  small  tension  in  the  link  adjacent  to  an  inserted  node.  This  abrupt 
error  then  destabilized  the  computation.  The  variable  length  cases  were  revised 
so  that  no  nodes  were  inserted  during  a  run.  Other  users  have  not  reoorted  this 
problem  for  "full  size"  ocean  cable  systems. 

The  second  premise  for  SNAPLOAD  is  that  it  trades  computation  economy  for 
reasonable  accuracy.  However,  the  algorithm  for  controlling  the  time  step  lacks  a 
minimum  step  size.  There  is  the  risk  in  any  run  that  an  instability  early  in  the 
run  will  stall  the  computation  until  the  computer  time  allotted  for  the  entire 
case  is  exhausted.  This  can  be  a  substantial  waste  for  long  events. 

Although  SNAPLOAD  obtained  solutions  that  modeled  the  essential  events  with 
creditable  accuracy  in  most  cases,  its  potential  for  widespread  use  is  limited 
until  these  drawbacks  have  been  addressed. 

One  of  the  attractions  for  SNAPLOAD  users  is  that  it  is  more  economical  than 
SEADYN.  Even  in  the  variable  length  cases  reported  here,  where  SNAPLOAD  was  re¬ 
quired  to  converge  to  extremely  narrow  error  bounds,  it  used  only  one-fifth  to 
one-third  the  computer  time  required  by  SEADYN  for  equivalent  runs.  One  may 
expect  a  greater  cost  advantage  for  SNAPLOAD  for  "easier"  cases. 
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RECOMMENDATIONS 


This  series  of  four  experiments  provides  a  basis  for  evaluating  dynamic 
models  of  ocean  cable  structures.  The  data  have  been  checked  by  comparison  to 
elementary  theory  at  several  points.  The  data  base  includes  sizes  of  6  feet, 

60  feet,  as  well  as  variable  length  and  the  deployment  of  a  scientific  mooring  in 
the  ocean.  The  data  base  should  be  extended  to  include  the  massive  components 
used  in  ship  moorings  and  other  large  ocean  structures.  The  value  of  dynamic 
modeling  with  computers  is  not  realized  for  lightweight  moorings.  Some  specia¬ 
lized  scientific  moorings  may  require  sophisticated  "tuned"  iiioorings  in  order  to 
make  reliable  measurements,  but  the  major  contribution  of  dynamic  modeling  is  to 
reduce  the  risk  of  loss  or  damage  to  enormously  expensive  large  ocean  structures. 
It  is  essential  that  dynamic  models  be  verified  for  full-size  ocean  structures. 

The  models  studied  in  this  report  have  shown  that  modeling  the  elastic  pro¬ 
perties  of  cables  has  an  important  effect  on  the  stability  of  the  models,  even  in 
cases  where  the  elastic  response  of  the  cable  has  only  a  slight  effect  on  the 
dynamic  event  itself.  Load  versus  strain  data  for  a  wide  variety  of  cable  con¬ 
structions  may  be  available  at  modest  cost  from  manufacturers  willing  to  release 
sample  results  of  their  routine  quality  assurance  load  testing.  Material  damping 
data  are  scarce,  and  a  program  to  measure  the  damping  of  various  materials  and 
cable  constructions  should  be  developed  and  implemented.  The  concept  of  the 
variable  length  experiment  is  suggested  as  a  technique  for  finding  both  the  effec¬ 
tive  elasticity  (EA)  and  material  damping  factor.  The  limited  experience  gained 
in  these  few  cases  is  impressive  because  the  results  correlated  well  with  elemen¬ 
tary  spring-mass  theory  even  though  the  materials  used  were  nonlinear  plastics. 
Furthermore,  the  method  allows  water  absorbing  materials  like  nylon  to  be  immersed 
during  the  test  measurement. 

Both  models  have  shown  that  they  accurately  compute  the  major  features  of  an 
event.  For  both  models,  however,  there  is  a  clear  need  for  a  continuing  program 
of  maintenance.  This  does  not  require  full-time  attention,  but  it  is  essential 
that  a  protocol  for  correcting,  completing,  and  adapting  these  codes  be  main¬ 
tained.  Both  programs  have  minor  errors  that  should  be  corrected.  SNAPLOAD  has 
lacked  adequate  maintenance  in  the  past. 

The  time  step  algorithm  in  SEADYN  needs  to  be  adapted  to  recognize  the  onset 
of  a  "mitosis"  event  during  variable  length  modeling. 
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The  SEADYN  manual  is  oriented  towards  the  parameters  that  the  user  inputs  to 
define  a  case  and  the  way  SEADYN  will  process  it.  Each  parameter  is  listed  and 
described  in  the  manual.  These  descriptions  should  be  expanded,  giving  the  user 
more  guidance  in  the  selection  of  values  to  use.  This  may  rely  on  cataloging  what 
value  has  worked  with  which  prior  case. 

An  archive  of  SEADYN  input  decks  for  successful  problems  should  be  kept,  along 
with  annotations  of  critical  parameters.  A  copy  of  this  archive  should  accompany 
the  manual  sent  to  new  users. 

SEADYN  is  a  very  general  model,  both  in  its  mathematical  formulation  and  in 
the  set  of  options  programmed  for  it.  This  presents  a  formidable  obstacle  to  the 
novice  user.  Once  the  basic  algorithms  have  been  verified  for  a  suitable  variety 
of  problems,  it  is  recommended  that  a  series  of  programs  be  written  that  apply  the 
verified  method  to  smaller  classes  of  problems.  One  cf  these  should  be  serially 
connected  cables  in  a  vertical  plane,  because  this  is  a  very  common  ocean  struc¬ 
ture.  There  are  several  advantages  to  having  a  suite  of  smaller  programs  to 
augment  the  general  code.  First,  the  directed  purpose  reduces  the  number  of 
options  in  the  code.  This  results  in  simpler  input,  which  means  the  manuals  are 
easier  to  understand  (as  well  as  write).  The  code  itself  is  more  efficient,  not 
only  because  it  may  be  restricted  to  planar  structures,  but  also  because  the  time 
step  and  convergence  limits  can  be  optimized  for  the  kind  of  problem  addressed. 

There  are  several  recommendations  about  drag  coefficients: 

•  Replace  the  present  in-line  default  drag  coefficient  functions  with 
calls  to  subroutine  DRAGCO; 

•  Select  new  default  functions  that  more  closely  predict  the  drag 
forces; 

•  Rewrite  subroutine  DRAGCO  in  a  menu  format.  That  is,  provide  a  list 
of  drag  coefficient  versus  Reynolds  Number  functions  (starting  with  A 
+  B  *  (R**C)).  Calling  arguments  include  the  function  number  as  well 
as  the  function  coefficients.  The  above  function  will  handle  the  most 
common  cases  (A  +  B/R,  and  B  *  R**C)  by  setting  C  *  -1  or  A  =  0, 
respectively.  The  default  function  number  and  its  coefficients 
should  be  programmed  in  DATA  statements  using  dummy  arguments  to 
prevent  erasure. 

•  Provide  additional  input  cards  for  each  cable  material.  Each  card 
provides  the  drag  coefficient  function  number,  function  coefficients 
A,  B,  C,  etc,  as  well  as  the  Reynolds  Number  range  limitation.  A 
field  on  the  cable  material  card  itself  would  show  how  many  drag 
function  cards  to  read.  A  zero  entry  in  this  field  would  indicate  the 
default  functions.  Separate  cards  would  be  needed  for  normal  and 
tangential  functions. 
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•  Some  cables  generate  a  third  component  of  hydrodynamic  force,  trans¬ 
verse  to  the  cable  plane.  Addition  of  this  detail  would  be  convenient 
during  this  revision  of  the  code. 

The  foremost  recommendation  for  the  SNAPLOAD  program  is  that  a  protocol  be 
established  to  maintain  and  update  the  code.  Initial  tasks  include  the  following: 

•  Add  an  interpreter  in  the  input  and  output  stages  sc  that  nodes  count 
by  1  and  the  ambiguity  between  links  and  segements  is  resolved. 

•  Provide  a  minimum  time  step  or  other  "trigger"  to  stop  the  program 
when  the  algorithm  gets  stuck. 

•  Remove  the  170  output  table  limit  by  using  an  extra  I/O  unit, 

•  Adjust  the  code  so  that  payout  nodes  are  inserted  smoothly,  even  for 
laboratory-scale  cases. 

•  Add  material  damping  to  the  code  or  modify  the  mathematical  damping  so 
that  it  relates  to  the  physical  damping  property. 

•  Add  a  new  algorithm  to  calculate  the  initial  cable  equilibrium  state 
accurately  from  an  initial  user  definition  as  a  straight  horizontal  or 
vertical  line. 

•  Change  the  default  drag  coefficients  and  add  a  drag  coefficient  sub¬ 
routine  as  recommended  for  SEAOYN. 

•  Make  other  trivial  parameters  that  are  fixed  constants  in  the  present 
code  user  input  variables.  The  fluid  density  is  an  example;  SNAPLOAD 
models  cables  suspended  in  seawater  only. 

Finally,  the  material  damping  parameter  has  had  such  a  striking  affect  in 
smoothing  the  tension  response  for  the  variable  length  comparisons,  that  certain  of 
the  comparison  cases  reported  here  should  be  remodeled.  In  particular,  the  dual 
cable  suspended  load  run  133  from  the  6-foot  series,  the  deployment  simulation 
run  6  from  the  60-foot  series,  and  the  MDE  deployment  were  badly  affected  by  the 
lack  of  material  damping  in  the  initial  comparison.  The  remodeling  will  verify 
this  assertion.  In  addition,  remodeling  the  6-foot  case  will  verify  that  the 
development  of  SEADYN  since  those  comparisons  were  performed  has  been  consistent 
with  the  earlier  results. 
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A . 1  ASSIGNING  THE  NODES. 

Given  an  overall  length  of  cable  in  a  system,  the  user  must  select  the  number 
and  location  of  nodes  to  be  included  in  his  input  model.  There  will  usually  be  a 
fairly  obvious  minimum  number  of  nodes,  located  at  cable  ends  and  Y's  and  where 
discrete  bodies  are  attached.  But  more  nodes  will  usually  be  required  in  order  to 
express  the  dynamic  curvature  of  the  cable  adequately. 

However,  there  is  a  double  penalty  extracted  for  using  more  nodes.  On  the 
one  hand,  more  nodes  mean  more  computation  in  each  "pass"  along  the  cable.  On  the 
other  hand,  more  nodes  mean  shorter  elements.  There  is  a  direct  relationship 
between  element  length  and  time  step  size.  Shorter  elements  require  shorter  time 
steps.  The  number  of  time  steps  required  to  model  a  given  time  interval  therefore 
increases  with  the  number  of  nodes  used  in  the  model.  When  the  user  defines  more 
nodes,  the  computer  executes  more  passes  requiring  more  computation. 

The  user  relies  on  experience,  intuition,  and  the  results  of  trial  runs  to 
concentrate  nodes  in  areas  of  sharpest  curvature. 

Figure  2-21  shows  the  CEL  mooring  as  it  was  deployed  during  the  MDE.  Fig¬ 
ure  4-46  shows  the  mooring  as  it  was  modeled  for  SEADYN.  By  comparing  these 
figures,  it  will  be  seen  that  nodes  were  located  where  one  or  more  instruments 
were  clustered.  Long  wire  rope  spans  were  broken  into  uniform  elements  between 
100  and  200  feet  long. 

The  center  spans  were  assigned  longer  elements  then  the  spans  near  the  upper 
and  lower  ends,  because  the  cable  curvature  is  greatest  near  the  ends  in  this 
problem.  Element  lengths  less  than  100  feet  were  accepted  when  they  could  not  be 
avoided,  namely,  elements  1,  2,  and  21  where  the  node  spacing  defies  the  element. 
Element  lengths  were  made  long  enough  that  the  transit  time  for  a  tension  wave 
along  an  element  is  at  least  0.01  seconds.  This  allows  a  0.005  second  time  step 
to  be  used  in  SEADYN  without  numerical  instability. 

A. 2  NODAL  DRAG  COEFFICIENT. 

The  terminal  velocity  of  the  combination  of  nodes  1,  2,  and  3  is  shown  on 
Figure  A-1  by  the  slope  of  the  depth  trace  for  pinger  ALPHA.  During  the  first 
2  minutes,  the  descent  speed  is  6.14  feet  per  second.  Later  in  the  descent,  the 
drag  and  buoyancy  of  other  nodes  slow  the  descent  perceptibly. 
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FIGURE  A-1  MDE  Experiment  5  PINGER  A  Burst  2  Depth 


From  Table  A-1,  the  net  weight  of  nodes  1,  2,  and  3  is  2600-611-79C  =  1199 
pounds.  Using  a  familiar  hydrodynamic  expression  in  the  form 

Cp.A^  =  F/Q, 

where  F  is  force  and  Q  is  dynamic  pressure, 

gives  ^D‘^T  ” 

The  actual  anchor  was  made  of  sandbags  heaped  Oii  a  wooden  pallet  and  con¬ 
tained  by  a  coarse  net.  This  irregular  shape  was  approximated  by  a  sphere  4  feet 
in  diameter.  The  diameter  of  the  sphere  was  selected  to  coincide  with  the  length 
of  a  side  of  the  pallet.  The  spherical  shape  was  assumed  to  approximate  the  shape 
of  the  sandbag  pile.  SEADYN  is  able  to  calculate  the  added  mass  coefficients  for 
spheres.  Combining  the  frontal  area  of  the  anchor,  assumed  to  be  a  4-foot  dia¬ 
meter  sphere,  with  the  areas  of  the  components  in  nodes  two  and  three  gives 

Ay  =  12.4  +  16.3  +  12.1  =  41  ft^, 

so  that  Cg  =  32/41  =  0.78. 

It  is  assumed  that  this  coefficient  applies  to  all  the  nodes,  not  only 
nodes  1,  2,  and  3  as  computed. 

A. 3  NODAL  WEIGHT. 

The  weight  of  each  node  is  simply  the  sum  of  the  weights  of  the  parts.  The 
weights  and  dimensions  of  the  MDE  components  were  taken  from  Table  2-2.  Table  A-1 
shows  how  these  weights  were  assigned  to  each  node. 

An  element  in  SEADYN  extends  from  the  center  of  one  node  to  the  center  of  the 
next.  As  shown  on  Figure  2-16,  the  actual  instruments  were  inserted  in  the  MDE 
mooring  line,  not  clamped  alongside.  Thus,  the  element  lengths  tabulated  in 
Figure  4-46  overlap  the  length  of  the  instruments  themselves.  The  weights  of  the 
nodes  are  adjusted  downward  to  account  for  the  weight  of  this  extra  length. 

A. 4  EQUIVALENT  DIAMETERS. 

Table  A-1  also  includes  the  frontal  areas  of  each  nodal  component.  The 
equivalent  diameter  for  a  spherical  node  is  calculated  to  give  the  same  frontal 
area.  Note  that  this  diameter  must  not  be  used  to  calculate  the  "equivalent 
displacement"  volume.  That  is  why  the  weights  used  in  SEADYN  are  immersed  values. 
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TABLE  A-1.  MDE  NODAL  PARAMETERS 


NODE 

COMPONENT 

WEIGHT 

(LB.) 

FRONTAL  AREA 
(FT2) 

EQUIVALENT 

DIAMETER 

(FT) 

2 

Release  1 

70. 

2.14 

Release  2 

70. 

2.14 

Frame  (estimated) 

20. 

2.00 

OC  Pinger 

76. 

2.20 

28"  sphere 

-847. 

7.88 

. 

m. 

lOS" 

4.56 

3 

38"  Sphere 

-847. 

7.38 

DC  Pinger 

76. 

2.20 

FVR 

-4. 

1.97 

Chain  Excess* 

-15. 

-790. 

T05 

3.92 

8,  13 

SC  Pinger 

44.3 

1.15 

FVR 

-4. 

1.97 

Chain 

10.2 

.51 

Wire  Excess* 

-.5 

3.63 

2.15 

11 

T/P 

Wire  Excess* 

.96 

21 

SC  Pinger 

1.15 

FVR 

1.97 

T/P 

20.0 

0.72 

Wire  Excess* 

-.5 

5575” 

3.84 

2.21 

22 

44"  Sphere 

-1050. 

Chain  Deficit* 

4. 

T55F 

3.67 

*  Corrects  for  excess/deficit  in  elenent  lengths  adjusted  to  reach  center  of 
spherical  node. 
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INPUT  DATA  FOR  SEADYN  AND  SNAPLOAD 
VARIABLE  LENGTH  COMPARISONS 


APPENDIX  B 


INPUT  DATA  FOR  SEADYN  AND  SNAPLOAD 
VARIABLE  LENGTH  COf^PARISONS 


Tables  B-1  through  B-6  show  input  data  decks  for  the  SEADYN  comparison  runs. 
The  descriptions  of  the  card  images  correspond  to  major  card  description  headings 
in  the  SEADYN  user's  manual  (Reference  10).  Tables  B-7  through  B-9  show  FORTRAN 
update  decks  for  the  SEADYN  source  code  subroutine  TVARY.  This  subroutine  de¬ 
scribes  the  time  varying  moving  boundary  and  payout  velocities. 

Tables  B-IO  through  B-15  show  input  data  decks  for  SNAPLOAD.  Reference  11 
contains  the  SNAPLOAD  user's  guide. 

Table  B-16  is  an  input  deck  for  SEADYN.  It  is  used  to  estimate  the  material 
damping  parameter.  The  event  models  a  cable  slightly  stretched  between  fixed  sup¬ 
ports.  The  center  node  is  displaced  slightly  and  abruptly  released  at  the  start. 
The  material  damping  parameter  is  adjusted  in  successive  runs  until  the  oscilla¬ 
tion  of  the  center  node  decays  to  half  its  initial  displacement  in  a  certain 
number  of  cycles. 
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TABLE  B-2.  SEADYN  INPUT  DECK  FOR  VARIABLE  LENGTH  TEST  12 
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TABLE  B-4.  SEADYN  INPUT  DECK  FOR  VARIABLE  LENGTH  TEST  29 
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TABLE  B-7.  SEADYN  SUBROUTINE  TVARY  STATEMENTS  FOR 
VARIABLE  LENGTH  TESTS  3  AND  12 


*IDENT  LK1024 

C  TUARY  UPDATES  FOR  TEST  3 
*11  TVARY.  6* 8 

IF(I.ECJ.2)  GO  TO  600 
IF<T.GE.O.O. AND.T.LE.2,2) 
IF<T.GT.2.2.AND.T.LE.23.2) 
IF ( T . GT . 23 . 2 . AND . T . LE . 25 . 4 
IF<T.GT.25.4)  F=0. 

RETURN 

600  F=COS( 15.708*T) 

RETURN 
100  CONTINUE 


*IDENT  LK1024 

C  TVARY  UPDATES  FOR  TEST  12 
*D  TVARY.6f8 

IF<I,EQ.2)  GO  TO  600 
F=0. 

IF<T.GT,4.0.AND.T.LE.6. )  F 
IF(T.GT.6.0.AND.T.LE.29.2) 
IF(T.GT.29.2.AND.T.LE,29,8 
IF(T.GT.29.8)  F=0. 

RETURN 

600  F=COS( 15.708»T) 

RETURN 

100  CONTINUE 
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=.8136*T 
F=1 .78992 
F=.8136*(25.4-T) 


=.8925*T-3.57 
F=1 .785 

)  F=2.98333»(29.8-T) 


TABLE  B-8.  SEAOYN  SUBROUTINE  TVARY  STATEMENT  FOR 
VARIABLE  LENGTH  TESTS  23  AND  29 


*IDENT  LK1024 

C  TVARY  UPDATES  FOR  TEST  23 
»D  TVARY.AfB 

IF(I.EQ.2)  GO  TO  600 
F=0. 

IF(T.GT.4.0.ANli.T.LE.6.  >  F= . 8925*T-3 , 57 
IF<T.GT.6.0.ANli.T.LE.30.4)  F  =  1 .785 
IF(T.GT.30.4.AND.T.LE.31 .2)  F=2 . 23125* ( 31 . 2-T ) 
IF<T.GT.31 .2)  F=0. 

RETURN 

600  F=COS< 15.708*T) 

RETURN 

100  CONTINUE 


*IDENT  LK1024 

C  TVARY  UPDATES  FOR  TEST  29 
*D  TVARY. 6»8 

IF<I.EQ.2)  GO  TO  600 

IF(T.GE.O.O.AND.T.LE.  .8)  F=2.5''*T 

IFCT.GT. .8.AND.T.LE. 13.5)  F=2.056 

IFCT.GT. 13.5.AND.T.LE.14.0)  F=4 . 1 12* ( 1 4 . 0-T ) 

IF(T.GT.14.0)  F=0. 

RETURN 

600  F=COS( 15.708*T) 

RETURN 

100  CONTINUE 


table  B-9.  seadyn  subroutine  tvary  statement 
FOR  variable  length  TESTS  47  AND  48 


♦IDENT  LK1024 

C  TVARY  UPDATES  FOR  TEST  47 
*D  TVARY.6.a 

IF<I.EQ.2)  GO  TO  600 
F=0. 

IF<T.GT.440.AND.T.LE.6. >  F=T-4. 

IF(T. GT, 6.0. AND. T.LE. 17.6)  F=2. 
IF<T.GT.17.6.AND.T.LE.18.1)  F=4  .» ( 18 . 1-T ) 
IF<T.GT.18.1)  F=0. 

RETURN 

600  F=C0S<9.4248*T) 

RETURN 
100  CONTINUE 


«IDENT  LK1024 

C  TVARY  UPDATES  FOR  TEST  48 
*D  TVARY.6»8 

IF<I.EQ.2)  GO  TO  600 
F=0. 

IF(T.GT.4.0.ANEi.T.LE.6.)  F“T'”4. 
IF(T.GT.6.0.AND.T.LE.18.2)  F=2. 

IF (T.GT, 18. 2. AND. T.LE. 18. 7)  F=4.#( 18. 7-T) 
IF<T.GT. 18.7)  F=0. 

RETURN 

600  F=-SIN(9.4248*T> 

RETURN 
100  CONTINUE 
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TABLE  B-11.  SNAPLOAD  INPUT  DECK  FOR  VARIABLE  LENGTH  TEST  12 
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TABLE  B-12.  SNAPLOAO  INPUT  DECK  FOR  VARIABLE  LENGTH  TEST  23 
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TABLE  B-14.  SNAPLOAD  INPUT  DECK  FOR  VARIABLE  LENGTH  TEST  47 
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TABLE  B-16.  SEADYN  INPUT  DECK  FOR  ESTIMATING  MATERIAL  DAMPING 


o 

C-J 


O  UT  •H  o  ^ 
'O 
fs 
o 

^  ^  o 


\n  yn 


C3 

S3 

t-.  ^  o  ^  n  C  l 
O  I  I 

o 


o 

o 


o 

o 


o 

GO 

rs. 

ro 

ri 

fv  o 
»• 
CD 
N 

o  y) 

ri  tfi 

o  ^ 

•  n 

Cl 

y3  o 

GO 

y) 

U") 

«r 

ro 

Cl 

\n  o 
> 
03 

IS 

'O 

\n 

«T 

ro 

Cl 

^  o 
> 

fs 
'O 
jTi 

n 

ri 


n 


z  o 

3 
iZ 

a 


C4 


c 

O  ^  I 


O  O  S'  O  O  o  W 
o  o  o  c 

O  O  O  Q 

•  •  ^  o 

I  o  ^  « 

^  *-4  S'  ^ 

Cl 

9* 


!»■) 
o  o 
o  o 
o  \n 

o  • 

♦  Cl 

\n  Cl 


Cl  n 

fS  1-4 

o  ^  o  ^  v-t 

\n 

o 

o  rs  o  o 

o 

• 

o  o  o  o  U1 

o 

fl  •  •  • 

•  o  •  o  o 

o  • 

n  o  o  o 

•  o 

o  o 

^  Cl 

o  ^  «-4  no 

oc  M  f  I  Cl  n 

UJ 


< 

z 


^  Cl  ri 

o  o  *>4  a  1*4 
O  C  O  4C  o 
•  Ul  O  UJ  o 


o 

^  llT 

z  o  o  a 
>-  o  o  z 
a  •  •  ui 


B-17 


